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ABSTRACT 

We present the metalhcity distribution function (MDF) for 24,270 G and 16,847 K 
dwarfs at distances from 0.2 to 2.3 kpc from the Galactic plane, based on spectroscopy 
from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. 
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This stellar sample is significantly larger in both number and volume than previous 
spectroscopic analyses, which were limited to the solar vicinity, making it ideal for 
comparison with local volume-limited samples and Galactic models. For the first time, 
we have corrected the MDF for the various observational biases introduced by the 
SEGUE target selection strategy. The SEGUE sample is particularly notable for K 
dwarfs, which are too faint to examine spectroscopically far from the solar neighborhood. 
The MDF of both spectral types becomes more metal-poor with increasing |Z|, which 
reflects the transition from a sample with small [a/Fe] values at small heights to one 
with enhanced [a/Fe] above 1 kpc. Comparison of our SEGUE distributions to those 
of two different Milky Way models reveals that both are more metal-rich than our 
observed distributions at all heights above the plane. Our unbiased observations of G 
and K dwarfs provide valuable constraints over the |Z|-height range of the Milky Way 
disk for chemical and dynamical Galaxy evolution models, previously only calibrated 
to the solar neighborhood, with particular utility for thin- and thick-disk formation 
models. 

Subject headings: astronomical databases: miscellaneous - astronomical databases: sur- 
veys - Galaxy: abundances - Galaxy: disk - Galaxy: evolution - Galaxy: formation - 
Galaxy: stellar content - Galaxy: structure - stars: abundances - stars: distances 

1. Introduction 

Measuring the metallicity distribution of stars in the Milky Way disk is imperative for under- 
standing its chemical and dynamical evolution. Cool stars, such as G and K dwarfs, have lifetimes 
comparable to the age of the Galaxy, providing a complete fossil record of chemical development. 
By measuring the metallicity distribution of these cool stars, we provide constraints on the disk's 
star-formation history, and how it varies with respect to time and location. The metal-poor end 
of the metallicity distribution function (MDF) reveals information about the earliest era of star 
formation in the Galaxy, such as the mass function of the earliest stars, the extent of chemical 
pre-enrichment, and the rates and yields of core collapse supernovae. The metal-rich end of the 
MDF reflects recent Galaxy conditions, such as the present-day mass function and the frequency 
and yields of both type la and II supernovae. The MDF provides information about the regulation 
of star formation and the merger and accretion history of the Galaxy. 

There is a wide range of chemical and dynamical evolution simulations that model the structure 
of the Galactic disk. These vary with respect to properties such as the star-formation history, 
initial mass function, prominence of inflows and outflows, and the likelihood of mergers. All use 
observed samples to test their predictions of the metallicity distribution. Until recently, we lacked an 
adequate observational sample to quantitatively test these models beyond the solar neighborhood. 
In this work, we use the Sloan Extension for Galactic Understanding and Exploration (SEGUE, 
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Yannv et al.ll2009l ) survey to determine an unbiased MDF of cool dwarfs over a large volume of the 
disk, from 0.2 to 2.3 kpc from the plane of the Galaxy. 



I nitial examination of the c hemistry of local G dwa r fs revealed the "G-dwarf probl e m" (Ivan den Bergh 



19621 : IPagel fc Patchettl Il975l : IWvse fc Gilmorel Il995l : iRocha-Pinto fc Maciei Il99d : IPavata et al 



19971). Earl y mod els of star formation and chemical evolution, such as the simple closed box model 



of ISchmidtl (jl963l ). predicted many more low-metallicity G dwarfs than were actually observed. 
It was initially suspected that the G-dwarf problem arose from observational biases, namely that 
metal-rich stars are brighter and thus likely to be over-represented in a magnitude-limited sam- 
ple. However, the deficiency of low-metallicity stars persisted in later observational samples that 
corrected for these biases. Work such as the Geneva-Copenhagen survey ( GCS) of the solar neigh 



2011 



borhood (jjorgenseni l200d : iNordstrom et al. 



2004 



Holmberg et al. 



2007 



2009 



Casagrande et al. 



) indicated that for a large volume-complete and kinematically-unbiased sample, the simple 
closed box model over-predicted the number of cool metal-poor dwarfs even more than originally 
reported. 

We know that many of the assumptions in the simple closed box model are unphysical, in 
particular instantaneous recycling and the absence of gas flows. More recent models have updated 
and experimented with the model parameters to better mat ch the observed metallicity distribu - 



tion, for example, adding inflows of low-metallicity material (jLarson 



varying the formation time scales of different Galaxy components (IChiappini et al. 



mencing star formation from pre-enriched material (ITruran &: Cameron 



experimented with the initial mass function (IMF) (jChiappini et al 



1972: 



Chiappini et al.l 119971 ) 



2001 



2001 



or com- 



19711). Models h ave also 
; iRomano et al.l l2005l ) ; a 



proposed metallicity-dependent IMF will produce more high- than low- mass stars at early times, 
resulting in fewer cool metal-poor stars. 

This exploration of the G-dwarf problem prompted observations of the MDF of cooler stars. K 
and M dwarfs have even longer lifetimes than G-dwarf stars. There is a chance that more metal-rich 
G dwarfs have evolved off of the main sequence, leading to a metallicity bias in the MDF, which will 
not occur for cooler spectral types. In addition, comparing the MDFs of different spectral types is 
particularly useful for constraining the variation in the IMF. For example, if more low-mass stars are 
created as the metallicity of the enviro nment increas e s, there should be re l atively fewer metal-poo r 
K dwarfs than G dwarfs. Work s uch aslMouldl (Il982l) . IPavata et all (Il997l) . iFlvnn k M(^ (Il997l ). 



Rocha-Pinto k Maciell (|l998l ). iKotoneva et al.l (120021 ). and IWoolf k WestI (120121 ). found that, just 
as with G-dwarf stars, the simple closed box model predicted more metal-poor stars than were 
observed for both K and M dwarfs, implying that the IMF and star-formation history for these 
different spectral types are similar to one another. 

Measurements of the G- and K-dwarf MDF have been mostly confined to the solar neighbor- 
hood, limiting our understanding of disk properties and evolution processes to the local volume. 
However, the chemical and dynamical structure of the disk does not appear to be uniform with 
respect to Galactic height, but rather composed of two components, a thin and thick disk, of 
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unclear interdependence. iGilmore &: ReidI (|l983l ) first detected the thick disk in the Milky Way 
when they determined that the stellar number density as a function of height above the plane was 
best fit by two component^, one with a scale height of approximately 300 pc, and the s econd 
with a scale he i ght of 1350 pc. Analyses bv lGilmore &: Wvsd (|l985l ): IWvse &: Gilmord (jl995l ). and 
Chiba k BeeT] (|200ol ) established that the two populations were distinct in [Fe/H]. The so-called 
thick disk was metal-poor, with a peak metallicity around [Fe/H]= —0.6, in contrast to the thin 
disk, with a peak metallicity typically around [Fe/H]= —0.2. Further observ ation of solar neighbor - 
hood samples revealed that the two popula tions were distinc t in kinematics (ISoubiran et al.ll2003l) 
age ( Fuhrmann 19981). and a -abundance ( Fuhrmann 19981 : Prochaska et al. 200C : Bensbv et al 



2003, 



2OO5I : iReddv et ahlbood). Howev er, iNorris fc RvanI (|l99ll ) and recent work by Ibovv et"al 



(|2011al lbl) and iLiu Sz van de VenI (j2012l ). focusing on [a/Fe] vs. [Fe/H] and kinematics, questioned 
whether or not the two components were actually separable from one another, proposing that 
the Milky Way disk has a "thicker disk component," rather than two distinct structures. The 
debate over the basic structure of the Milky Way disk emphasizes the need for an unbiased spec- 
troscopic sample that extends over a large volume of the disk, such that we can investigate the 
metallicity structure of both of the proposed components with a uniform large data set. Further- 
more, it is critical to use our own Galaxy to disentangle the mechanisms behin d disk develop- 
ment because the thick disk appears to be a regular featur e in numerous galaxies (|Bursteinlll979l : 



Dalcanton BernsteinI I2OO2I: Ivoachim fc DalcantonI [2008allbh . including those at redshifts as high 



as z~3 ([Elmegreen &: ElmegreenI 120061 ). In this work, we investigate and constrain the chemical 
structure of the Milky Way as a whole, including both thin- and thick-disk components. 



1.1. Previous Analyses of the Chemical Structure of the Disk 



Past observations of the MDF have been limited in accuracy, samp le size, and, most impor - 
tantly, volume. Due to t he low luminosity o f cool s tars, analyses such aslPaeel fc PatchettI (119751). 
Gilmore fc Wvsel (Il98,5l'l. IWvse fc Gilmorel (Il99,5l 'l Ifiocha-Pinto fc Maciei (Il996l). Ipivnu fc MorelJ 



Kotoneva et al 



(12002) relied on pho- 



(|l997l ). lR,ocha-Pinto fc Maciei (|l998l ). lje)rgensen[(|20od ). and 
tometric calibrations to determine metallicities. Photometric metallicity determinations are suscep- 
tible to errors from reddening corrections, have reduced sensitivity at low metallicity, and depend 
strongly on the adopted calibration to spectroscopic estimates, which vary from work to work. Us- 
ing spectroscopic measurements increases the accuracy and precision of metallicity determinations, 
in addition to providing kinematic information such as radial velocities, albeit with the significant 
added cost of increased observing time. Previous spectroscopic analyses of these low luminosity 



^T he likely existence of a second stellar component close to the Galactic plane was previously noted by lYoshii 
( 19821 ). wh o refer r ed to it as a halo component, even though its inferred density was 10 times that of the local 
halo. The IVoshiil |l982h normalization relative t o the local thin disk (0 .01-0.02) and scale height (~2 kpc) were 
commensura t e wit h the values later determined by Gilmore fc ReidI (1983) (0.02, and 1.5 kpc, respectively). See also 
Yoshii et all (| 19871 '). 
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targets were limited to hundreds of stars along individual lines of sight (IWvse &: Gilmorg 



Favata et al 



2007 



1997 



Fuhrmann 



Fuhrmann 



200. 



1998 



2004 



Arnadottir et al 



Allende Prieto et al 



20091 : 



Fuhrmann 



2011 



2004 



Luck k. Heitei 



Katz et al. 



2011 



2005 



1995 : 



2009 . 



Two recent surveys have improved analyses of the MDF of the Milky Way disk by increasing the 
sample size and sky coverage. The GCS is a magnitude-limi ted survey, with ^14,000 F and G dwarfs 



with metallicities estiniated from Stromgren photom etry (jJorgensenl l2000l : iNordstrom et al.l l2004l : 
Holmberg etaP l2007l . bood : ICasagrande et allboilh . The RAdial Velocity Experiment (RAVE) 
further expai ids the sample of c ool stars, with around 17,000 F and G dwarfs and spectroscopic 



metallicities (jSiebert et al.ll201ll ). Although they have uniform data sets over a large region of the 



sky, neither of these have dwar f stars far beyond the p lane of the Galaxy. The GCS sample is 

). The dwarfs in the RAVE sample probe 



limited to 200 pc from the Sun (ICasagrande et al 



2011 



to a maximum distance of approximately 1 kpc: the typica l distance for their cool dwarf sample 
ranges from 50-250 pc ([Zwitter et al.ll2010l : ISteinmeta |2012| ) . Previous analyses of individual line 
of sights were similarly limited, extending no further than 50 pc from the Sun. Both the GCS 
and RAVE samples exhibit M DFs that peak at aroun d solar metallicity, desp ite their different 

indicating that 



methods for estimating [Fe/H] (jCasagrande et al. 



2011 



Cogkunoglu et al. 



2011 



they are dominated by thin-disk stars, as expected for a solar neighborhood sample. Without a 
large sample of cool stars far beyond the solar neighborhood, we cannot accurately constrain the 
chemical and dynamical development of the Milky Way disk as a whole. 



In contrast to GCS and RAVE, the early work on G dwarfs of lGilmore &: Wysd (119851 ) reached 
Galactocentric height above the plane, |Z|, of around 1.6 kpc, beyond the local volume. This 
analysis, however, consisted of only ten lines of sight and estimated inetallic ity from ultraviolet 
excess, based on broadband photometry. More recent work bv lKatz et al.l (j201ll ) examines subgiants 
and giants as high as 5 kpc above the plane of the Galaxy. Although they use spectroscopic 
metallicities, this sample probes only ~400 stars along two lines of sight, limiting their statistical 
power to measure the extremes of the MDF or its variation over their sample volume. 

SDSS and SEGUE provide a uniform sample over a large portion of the Milky Way; previous 
analyses have sought to use this data to const r ain t he chemical structure beyond the solar neigh- 
borhood. Ilvezic et al.l (j2008l ) and I Bond et al.l (j2010l ) examined over 2 million F and G dwarfs in 
SDSS, utilizing broadband ugr photometry to estimate both the metallicity and distance for each 
star. This large sample probes far beyond the local volume, ranging from \Z\ of 0.5 to 7 kpc, and is 
dominated by thick-disk and halo stars. They find that the disk metallicity decreases from around 
[Fe/H] of —0.6 close to the plane of the Galaxy to plateau at around —0.8. They also d etect a clear 
break between the disk and halo stars, which have a mean metallicity of around —1.4 (iBond et al. 
2010l ). Although they have a large and unbiased stellar sample, uncertainties in the photometric 



metallicity calibration affect the absolute scale of the MDF and manifest as errors in the derived 
photometric parallax relationships. They also cannot exclude bi naries from the sample ; an unde- 
tected companion has a significant effect on stellar photometry (ISchlesinger et al.ll2010l ). Finally, 
although they analyze a range of stellar types, they focus their analysis on F stars, specifying a 
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color range of (g — r) from 0.2 to 0.4. F dwarfs have shorter hfetimes than G dwarfs, and are thus 
an incomplete sample of the disk chemical evolution. 

Avoiding redd ening and calibratioii uncer tainties, past studies utilized SDSS spectroscopy to 
measure the MDF. lAllende Prieto et al.l (j2006l ) examined a sample of 22,770 F and G dwarfs from 
SDSS Data Release 3, ranging from 1 to 8 kpc in \Z\. They focus their analysi s on more metal-poor 



stars, examining the metallicity differences between the thick disk and halo. lAllende Prieto et al 



(j2006l ) find that the MDF of G dwarfs between 1 and 3 kpc of the plane, associated with the thick 
disk, have a peak around [Fe/H]= —0.7, significantly more metal-rich than those above this height, 
which are halo stars and exhibit a peak [Fe/H] of —1.6. With the release of SEGUE in Data Release 
6, the sample size and coverage in Galactocentric height and radius have improved significantly 
since this work, allowing us t o more accurately determiii e the MDF of G dwarfs with respect to 
spatial position. Furthermore, lAllende Prieto et al.l ()2006l ) do not account for observational biases 
that originate from the SDSS target selection algorithm, which did not sample the stellar range of 
metallicity and age equally. 



The recent work of iLee et al.l (j2011bl ) utilize the SEGUE spectroscopic G-dwa rf sample to 



exami ne th e properties of the a-separated thin- and thick-disk components. Similar to llvezic et al 



(|2008l ). and lBond et al.l (|201Cl ). they find that the thick-disk component is more metal-poor than the 
thin-disk component, and dominates far from the plane of the Galaxy. Although they use up-to-date 
SSPP parameters, they do not examine the entire G-dwarf sample available in SEGUE, selecting 
only stars specific ally targeted by SEGUE as G dwarfs, rather than all stars that fulfill the criteria. 
Lee et al.l (l2011bl ) also do not account for the target selection biases in SEGUE, and are thus biased 
towards meta l-poor stars, although these observational biases are less significant in their sample 
than those of lAllende Prieto et al.l (j2006l ) . Most importantly, rather than examining the MDF of 
the disk as a whole, they examine the relative metallicity structure of two a-separated components. 
We examine the MDF of the SEGUE sample without making any chemical or kinematic separation 
of the populations. 



1.2. This Work 

With a large number of stars and extensive sky and volume coverage, the complete SEGUE 
survey provides an ideal sample to examine the chemical abundance distribution in the Galaxy 



(jYanny et al.ll2009l ). The SEGUE data consist of SDSS ugriz photometry and spectroscopy for 
240,000 stars over a range of 14< g <20.3 in ~3500 square degrees on the sky. There are 50,210 
and 26,834 SEGUE stars that fulfill the G- and K-dwarf photometric criteria, reaching Galactic 
distances from the plane, \Z\, of around 3.5 kpc and ranging from Galactic radius, R, of 5 to 13 kpc. 
Not only is this the largest spectroscopic sample available, it covers a much more extensive volume 
of the Milky Way disk than all previous analyses. Furthermore, it is the largest and most extensive 
spectroscopic sample of K-dwarf stars, whose faint magnitudes make it difficult to probe far from 
the solar neighborhood. Previous samples of K-dwarfs were limited in volume, spatial-coverage. 
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sample size, and atmospheric parameter accuracy. 

We utilize this sample to determine an unbiased MDF of both G- and K-dwarf stars throughout 
a large volume of the disk and determine how they change with respect to height above the plane. By 
using both spectral types, we examine a larger volume of the disk, constraining the disk chemical 
structure both close and far from the plane of the Galaxy, whereas analyses of hotter stars in 
SEGUE are limited to thick-disk dominated space. In addition, we compare the MDF of the two 
spectral types to analyze whether or not the star-formation and chemical-evolution history vary 
with respect to stellar mass. 

SEGUE selects stars as spectroscopic targets using a series of photometric and proper motion 
criteria to isolate stars of different spectral types. Because the SDSS and SEGUE survey addressed 
many different scientific questions, it had a large variety of target types. Each of these isolate a 
different portion of parameter space, resulting in samples which do not sample the stellar range of 
metallicity and age equally. The net effect for the entire sample is an observational bias in favor 
of metal-poor stars. This is the first work to systematically determine and account for the various 
observational biases in the SEGUE G- and K-dwarf sample. By comparing SDSS photometry and 
SEGUE spectroscopy for each line of sight, we assign weights to each SEGUE star to account for 
the various target-selection biases. Computing the MDF with this weighted sample allows us to 
examine how the underlying population of stars throughout a large volume of the disk varies with 
respect to spatial position. 

In this work, we first discuss how we extract the G and K dwarf stars from the SEGUE 
database (§[2]) and estimate the distance to each star using isochrones (§[3]). We then constrain the 
different observational biases that may arise in the SEGUE sample (§[3]) and detail our technique to 
correct for those biases that arise from the SEGUE target selection algorithm (§[121) ■ Utilizing these 
target selection weights, we determine complete and unbiased metallicity distribution functions for 
both G and K dwarfs that accurately reflect the underlying chemical structure (§[5]). In §[6l we 
compare the observed metallicity structure to previous analyses a nd two Galaxy models, the stellar 



p opulation synthesis rnodel T RILEGAL 1.4 (iGirardi et al.ll2005l ) and the radial migration model 



of lSch5nrich fc BinnevI I^OOgJH). These comparisons reveal how critical an accurate picture of the 
disk chemical structure is for understanding how the disk formed and evolved and how valuable 
our observations are for understanding how the Milky Way developed. 



2. The SEGUE Stellar Sample 

2.1. Basic Survey Design 

The SEGUE survey combines the extensive uniform data set of photometry from SDSS with 
medium-resolution (R~ 1800) spectro scopy over a broad spectral range (3800-9200A) for ~ 240,000 
stars over a range of spectral types (lYanny et al.ll2009l ). Technical information about the Sloan 
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Digital Sky Survey is published on the s urvey design dYork et al.1 boool: lEisenstein et al.l |2011| 



(Ivezic et al 



telescope and cam era (jGunn et alj 120061 . Il998l ) , astrometric (IPier et al.l 120031 ) and photometric 



2004) accuracy, photometr i c svstem (IFukugita et al.lll996l). and photoni etric calibra- 



tion (jHogg et al.ll200ll : ISmith et alJl2002l : lTucker et al.ll2006l : IPadmanabhan et al.ll2008l ). Beyond its 
large size, this sample has a homogeneous data set, substantial area coverage, and spectroscopically- 
detern iined stellar atmosphe ric parameters. This work utilizes photometry from Data Release 7 
(DR7. lAbazajian et al.ll2009l). The atmospheric parameters are a modified version of those released 
as part of DR8 (jAihara et al.ll201ll . Appendix lX|) . 

Each SEGUE plug-platei covers a circular region of 7 square degrees, probing the sky with 
640 spectroscopic fibers. SEGUE selects targets from SDSS for spectroscopic observation based on 
photometric and proper-motion cuts. For an individual SEGUE line of sight, approximately 375 
and 95 fibers are allotted to G and K dwarfs, respectively. The stars assigned spectroscopic fibers 
are randomly selected from all of the stars in SDSS photometry that meet the target-selection 
criteria. These criteria, and the number of fibers devoted to each target type, changed over the 
course of the SEGUE survey observations to improve the efficiency of some of the target selections. 
We use the target identifications from Data Release 7 applied uniformly to all of the photometric 
and spectroscopic data. We also limit ourselves to pointings which fall under the SEGUE program 
and not the extragalactic SDSS survey or the SEGUE-2 program, as neither of these targeted the G 
and K dwarfs explicitly. Finally, we eliminate any pointings that do not have both bright (rg <17.8) 
and faint (ro >17.8) plates; this requirement ensures we probe the same magnitude range for all 
lines of sight. 

The resulting spectra are processed through the SSPP, an automated system that determines 



atmospheric parameters, such as effective temperature, surface gravity, and metallicity (jLee et al 



2008al ) . The SSPP employs 6 primary methods for the estimation of Tefj, 10 for the estimation of 
log (7, and 12 for the estimation of [Fe/H]. For an in-depth description of all of the different SSPP 
calculations and techniques, see Lee et al. ( 2008al lbl). This program's outputs have been checked 



against high-reso l ution spectra of stars within globular and open clusters, as well as in the field 



(|Lee et al.ll2008bl : lAllende Prieto et al.ll2008l : ISmolinski et al.ll201ll ). The uncertainties of the SSPP 
for targets with S/N=25 per pixel, where each pixel is ~lA, are (T(Tcfr)=200 K, a{\og g)=0.4: 
dex, and (j([Fe/H])=0.3 dex. These uncertainties increase as the S/N de creases: for S'/A^=10, 



CT(Tcfr)=260 K, a{logg)=0.6 dex, and o-([Fe/H])=0.45 dex (ILee et al.ll2008al ) 



Not all of the techniques employed by the SSPP are accurate for analyzing cool dwarf stars. 
We have examined each method for estimating stellar parameters for our spectral types, refining 
and optimizing the SSPP for G and K dwarfs (see AppendixlAl). Specifically, for determinations 
of surface gravity and metallicity, we eliminate techniques which were designed for hotter and/or 
evolved stars and refine existing techniques to improve their accuracy. This revised SSPP was tested 
on open and globular clusters (see Table[l]). The modified techniques result in negligible shifts to 



^Throughout the rest of this work, we will refer to plug-plates simply as plates. 
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the overall metallicity determined by the SSPP for each cluster, well within the uncertainties. Our 
work on the revised SSPP will be part of the improvements to the SSPP released with SDSS DR9. 



2.2. Extracting G- and K-dwarf Stars from SEGUE 



SEGUE G and K dwarfs are selected using a simple color and magnitude cut. The SEGUE 
"G dwarfs" are defined as having 14. 0< tq <20. 2 and 0.48< (q - r )n <0.55, while the "K dwarfs" 
have 14. 5< ro <19.0 with 0.55< {g - r)o <0.75 (lYannv et al.ll2009ll4 Note th at the the subscript 



indicates dereddening and absorption correction using Schlegel et al. ( 19981) values. F or [Fe/H] 
from —0.5 to —2.5, the Yale Rotation Evolution Code isochrones (YREC. lAn et al.ll2009l ) indicate 
that these {g — r) colors correspond to a temperature range of ~4800— 5300 K for K dwarfs and 
~5000— 5600 K for G dwarfs. In this paper, we will refer to the SEGUE "G and K dwarf" categories 
simply as G and K dwarfs. These target categories were designed to isolate stars at a range of 
ages. Stars of hotter spectral types that formed early in the disk's history have evolved past the 
main-sequence turnoff, biasing samples towards younger, more metal-rich, stars. As we do not 
have parallax info rmation or precise surfa c e gravities, we can not e stima te the ages for our stella r 



sample, as done in lEdvardsson et al.l (jl993l ). lNordstrom et all (j2004l ) . and lCasagrande et al.l (j2011 
Fortunately, the effect of age on the MDF will be minimized for these cool stars. 

Extracting every star with SEGUE spectroscopy that matches these color and magnitude 
criteria from SDSS Data Release 7 photometry and the SSPP for Data Release 8 results in 50,210 
and 26,834 stars in the G and K categories, respectively. We then restrict our sample to targets with 
S/N >10, because thes e spectra have b etter-constrained and better-understood uncertainties than 
those with lower S/N (iLee et al.ll2008al ). We also eliminate targets where, for various reasons, the 
SSPP was unable to determine the temperature, metallicity, and/or surface gravity. In addition to 
catastrophic failures, we remove targets that the SSPP flags due to temperature or noise issues. For 
example, if the temperature determined for a star by the SSPP and that from a [g — z) relationship 
differ by more than 500 K, we eliminate it from our sample. Similarly, if a spectrum is flagged 
as noisy by the SSPP, we remove it from the sample, even if its reported S/N is greater than 
10. Finally, we select all stars with \ogg >4.1 to isolate dwarf stars. This surface-gravity cut is 
discussed in depth in ^ 14.31 We also adjust the magnitude limits of our sample. The SDSS saturation 
limit varies over the instrument. To ensure that the bright end of the sample is complete, we set 
a uniform bright limit of rg >15. In addition, for the sake of our observational bias corrections, 
we removed the faintest G-dwarf targets, with ro greater than 18.45 mag (see § l4.7.2p . With these 
parameter criteria, in conjunction with the color and updated magnitude limits, we have around 
26,600 G and 18,500 K dwarfs. The distribution of our G- and K-dwarf spectroscopic sample in 
various atmospheric parameters is shown in Figured} The values shown for each target are slightly 



^We fi nd that these c olor cuts are equivalent to a spectral type of G5 or G6 for G dwarfs and G7 through K2 for 
K dwarfs (Johnson 19631 ). 



different than those included in DR8 (jAihara et al.ll201ll ): specifically, the [Fe/H] and surface- gravity 
determinations have been optimized for the sample, as discussed in AppendixlAl 



3. Distance Determinations 



Distances are critical both for determining a volume-complete MDF and for examining the 
behavior of these stars over different regions of the Galaxy. As these stars lack trigonometric 
paralla xes, we est i mate distances by matching each star to YREC isochrones in color and metallicity 
(Fig.[2| lAn et alj 12003 ) . T here are numerous isoc hrone sets available, suc h as the Dartmout h 
(jDotter et alJbood ). BaSTI (jpietrinferni et alJbooj ). and Padova isochrones (joirardi et alJboO-d ). 
Each isochrone set will predict slightly different behavior along the main sequence, where our stars 
lie, resulting in systematic changes in distance between sets f^lB.4l) . We also estimate the distances 
using the photometric parallax relationship of llvezic et alJ (|2008l ) (S IB.Sp . 



We utilize the YREC isochrones for three main reasons. First, the YREC se t is directly 



calibrated on open and globular clusters, while the other isochrones are not. The llvezic et al 



(j2008) photometric parallax relationship is calibrated on the open and globular cluster distances 



from iHarrid (jl996l ) . Analysis in the interim has found significant var iation in these calculated 
distances, increasing the uncertainty of the calibration, e.g.. Table 2 in lAn et al.l (j2009l ). Second, 
as they are specifically designed to work with SEGUE observat ions, the YREC i sochrones are 
guided by clusters observed in the SDSS ugriz bands, whereas the llvezic et al.l ()2008l ) photometric 
parallax relationship must convert from Johnson-Cousins to ugriz, contributing additional distance 
uncertainty. Finally, the YREC isochrones cover the widest metallicity range of all the d i fferen t 



options, from [Fe/H] of —3.0 to +0.4. Both the other isochrones and the llvezic et al 
photometric parallax relationship extend to only [Fe/H] of —2.5. 



(120081 ) 



3.1. Isochrone Matching 

For each star, we select the closest YREC isochrone both above and below the target in 
[Fe/H]. We limit our isochrones to the main sequence, as we expect all of our targets to be in this 
evolutionary stage. We also assume an age of 10 Gyr; uncertainties associated with this value are 
discussed in ^ IB.5[ We then match the target's {g — r)o color to each of the bracketing isochrones, 
extracting the predicted parameters when the (^f — r)o color matches within 0.001 mag. Using 
linear interpolation, we determine relationships between stellar parameters, such as the absolute 
magnitude in ugriz, with respect to [Fe/H]. We then extract isochrone parameters using these 
relationships and the SSPP estimate of [Fe/H], applying the distance modulus to derive a distance 
in each filter. Our final value is the mean distance over all of the SDSS filters. We also estimate 
distances using cubic, rather than linear, interpolation, finding little difference between the stellar 
parameters determined from these two schemes. Figure[3] shows the distribution of distances for G 
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and K dwarfs. 



3.2. Distance Uncertainties 



There are a number of random and systematic uncertainties, such as the uncertainty in SSPP 
measurements and undetected binarity, that must be taken into account in our distance estimates. 
These uncertainties vary with respect to metahicity; smah changes in metallicity result in large 
changes in absolute magnitude at the met al-rich end, whereas there is a smaller effect at t he met al- 
poor end. We use the TRILEGAL 1.4 joirardi et al.l koO^ ) and ISch5nrich fc Binne^J j^OOgJH) 
Galaxy models (S I6.4.2I and l6.4.1]) to estimate the effect of errors in different parameters on our 
sample. These are discussed at length in AppendixiBl the total random and systematic distance 
errors are shown in Fig.[3l The random distance uncertainty originates in the photometric, SSPP 
[Fe/H], and [ce/Fe] errors. There is an additional random distance uncertainty that originates from 
isochrone choice. Although differences be tween each isoc hrone set will create a systemati c shift in 



distance, compariso n between th e YREC (lAn et al 



(jGirardi et al.ll2004l ) , and BaSTI (jCassisi et al.ll2006l ) isochrones reveals that each shifts the distance 



20091), Dartmouth (IDotter et al.ll2008l ). Padova 



in a different way. While distances determined with the Dartmouth set finds larger distances than 
YREC estimates, those from Padova and BaSTI are smaller. Thus, we treat this variation as a 
random uncertainty. We find that the total random distance uncertainty ranges from around 18% 
for stars with [Fe/H]> —0.5 to 8% for more metal-poor stars, dominated by uncertainty in the 
SSPP [Fe/H] determination. 

The systematic uncertainties in distance stem from our age assumptions and the possibility 
of undetected binarity. Both of these are discussed at length in AppendixiBl The most metal- 
rich stars are likely younger than our assumed age of 10 Gyr. This leads to a systematic shift in 
distance of —3% for the most metal-rich stars, while the metal-poor stars are largely unaffected. 
An undetected companion has a comparable effect but goes in the opposite direction; each pair 
will be systematically shifted by approximately +5% in distance over the en tire metallicity range. 
Cons ervatively, 65% of all G and K dwarfs are expected to have companions ([Puquennov &: Mavor 
Il99lh . leading to a total systematic shift of around +3% in distance due to undetected binarity. 

We adjust our distances to account for these systematic changes. We assign an offset to each 
star from age and binarity based on its optimized SSPP [Fe/H], which we then convolve with a 
Gaussian distribution. As our age assumptions will cause us to overestimate the distances, we 
subtract the estimated offset from our estimated distance to each star. Assuming each star is single 
will cause us to underestimate the distance to approximately 65% of the sample. We randomly 
select 65% of the sample and increase their distance by the assigned systematic offset from binarity. 
The systematic uncertainties are much smaller than the random uncertainties at the metal-rich end 
of the distribution, whereas, due to binarity, they are comparable at the metal-poor end. 
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3.3. Testing our Calculated Distances on Globular and Open Clusters 



None of our targets have measured parallaxes; thus, we cannot directly confirm that our 
calculated distances are correct. We instead test our methods against SDSS spectra of stars in 
open and globula r clusters, obser v ed for testing the SSP P determinations and SDSS photometry 
(jLee et al.l l2008bl : lAn et alJ l2009l : ISmolinski et alJ l201ll ) . Unfortunately, many of these studies 
focused on stars past the main-sequence turnoff because they are considerably brighter than dwarfs. 
We have targets along the main sequence for the clusters M13, M67, NGC 2420, and NGC 6791. 
This limits us to a metallicity range from [Fe/H] of —1.54 to +0.30. 



We select stars from the samples of iLee et al.l (j2008bl ) and ISmolinski et al.l (j201ll ) with well- 
determined [Fe/H], S/N >10, and logg >4.1 to ensure that our cluster members are on the main 
sequence. We do not cut on {g — r)o color for the globular clusters because this would severely 
limit our sample size. Furthermore, our distance calculations should hold for any star on the main 
sequence. As a check, we determined the average distance to the cluster derived from a color- 
constrained sample to that from the larger main-sequence sample; they match within the expected 
errors in distance. Thus, a color cut would not significantly affect our distance measurements to 
these clusters. 

T h e calc u lated clust e r dista n ces are listed in Table 21 in addition to other parame ters from 
HarrisI (jl996l ). Lee et al.l (|2008bl ). ISmohnski et al.l (|201ll ). and the WEBDA database (jPaunzen 
20081 ) . A comparison of the clusters and YREC isochrones is shown in FigurelH The observed 



stars agree w ell with the shifted YREC isochrones. The publi s hed distances to these clusters vary 
considerably (|Harrislll996l : iKraft &: Ivansll2003l : lAn et al.ll2009l : ISmolinski et al.ll201ll ): thus, we are 
pleased with the general agreement we observe from our isochrone-matching method. 



4. Constraining Possible Biases in the SEGUE Sample 

The selection criteria for G and K dwarfs create numerous biases which affect the sample. For 
our SEGUE G- and K-dwarf sample to accurately reflect the chemical properties of these spectral 
types in the Galactic disk, we must ensure that we constrain and account for these observational 
biases. In this section, we investigate and constrain the possible sources of contamination in the 
sample and the observational biases induced by the SEGUE target-selection algorithm. We also 
constrain the effect of each of these errors on our MDF. 



4.1. Photometric Errors 

Although photometric errors for SDSS photometry are small, typically 2-3% for each filter, 
these uncertainties may bump targets in and out of the G- and K-dwarf color and magnitude 
range. As targets become fainter, their photometric errors increase, making it more likely that 
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a fainter star will be shifted out of the SEGUE color and magnitude criteria due to photometric 
errors. 



Newbv et al.l (|201ll ) determine that the photometric uncertainties are constant up to rg of 
19.7. The SEGUE target-selection criteria for K dwarfs selects stars as faint as rg of 19, and are 
thus unaffected by the increasing photometric errors. The G-dwarf magnitude range extends to 
ro of 20.2; however, we trim our G-dwarf sample to rg < 18.45 in order to limit our sample to a 
magnitude range where the faint end is complete (^ I4.7.2p . The photometric uncertainties are thus 
consistent over the full magnitude range of our G- and K-dwarf sample. 

We examine the effect of photometric uncertainty on target selection using TRILEGAL models 
along SEGUE lines of si ght f^ | 6.4.2|]. We est imate the photometric errors in go and rg from the 



exponential functions of iNewby et al.l (120111 ) . We then convolve these errors with a Gaussian 
over multiple iterations to examine how the sample of G- and K-dwarf stars change with changes 
in photometry. Approximately 2% percent of stars are shifted in and out of the sample; those 
stars removed and inserted into the sample by photometric uncertainties cover a similar range in 
metallicity space. We examine how the MDF changes over the full metallicity range with this 
photometric uncertainty, comparing each iteration to the underlying MDF of the model itself. We 
find that the uncertainty in each bin of the MDF due to photometric errors is on average 6% and 
should not create any [Fe/H] biases in our derived MDF (see Fig. [5]). 



4.2. Undetected Binarity 



As the target-selection criteria for G and K dwarfs is based purely upon a color and magni- 
tude cut, quantifying t he effect of an undetected companion on photometry is critical. Work by 
Schlesinger et al.l ([20101) utilized numerical modeling of the SEGUE G- and K-dwarf stars to con- 
strain the effect of binarity on the photometry of the sample. A companion will typically change 
the (g — r) color measured for a target by 0.01ib0.02 magnitudes, less than the uncertainty in SDSS 
photometry. Undetected companions will also change the photometry in individual filters, affecting 
distance estimates (Appendix lB.6p . which we include in the total systematic distance uncertainty. 

Although this is a small change in color, it can still shift stars in and out of the {g — r) range 
for G and K dwarfs. Furthermore, as companions are typically cooler than the primary star, they 
will make targets appear redder, w hich may result in unde tected binarity preferentially removing 
the coolest stars from the sample. ISchlesinger et al.l (j2010l ) investigated the effect of binarity on 
the properties of SEGUE G and K dwarf s, finding that, assuming 65% of all G and K dwarfs 
are in binaries (jDuquennoy &: Mayorlll99ll ). the addition of a companion will bump around 1% of 
stars into and 2% of stars out of the {g — r) color criteria. Not only is this a minimal effect, but 
the amount of stars shifted in and out are roughly comparable, indicating that binarity will have 
little effect on the G- and K-dwarf target selection. Using TRILEGAL modeled lines of sight and 
a modeled population of companions (Appendix IB.Sh . we run a Monte Carlo analysis to estimate 
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the uncertainty in the MDF due to undetected companions. Binarity does not preferentiahy affect 
stars of a particular metahicity, with a typical uncertainty of 4% for each MDF bin. 



Beyond photometry, ISchlesinger et al.l ([20101) combine their numerical model with synthetic 



spectra processed through the SSPP to constrain how a secondary will affect the estimates of 
atmospheric parameters. Effective temperature, which we do not use in our analysis, is the most 
altered by undetected companions; around 18% of G- and K-dwarf primaries will be shifted by 
more than 150 K, the reported SSPP uncertainty for S/N of 50. However, undetected binarity 
has very little effect on the SSPP [Fe/H] estimates. Approximately 99% of the modeled G and K 
dwarfs undergo shifts of less than 0.2 dex in [Fe/H], less than the expected uncertainty in the SSPP 
estimates over all S/N. Thus, the effect of binarity on SSPP atmospheric parameters is well within 
the expected errors, and we do not take it into account. 



4.3. Subgiant Contamination 

The color and magnitude cuts used to identify our cool dwarfs do not exclusively target main- 
sequence stars. Subgiants and giants, specifically K giants, can fall into our sample and drastically 
affect the accuracy of our distance estimates, which assume all of the targets are on the main 
sequence. To isolate dwarf stars we apply a cut on the SSPP log 51, limiting it to 4.1 and above. 
However, the uncertainty of the SSPP surface gravity estimates for stars with S/N of 25 is ±0.4 
dex, which may result in evolved stars contaminating the dwarf sample. 

We use a three- fold approach to examine the extent of the subgiant contamination. First, we 
employ Galaxy models to estimate the number of subgiants expected to fall in our SEGUE G- 
and K-dwarf sample. Second, we manufacture a series of synthetic spectra at various evolutionary 
stages and process them through the SSPP. This allows us to determine if there is a particular 
region of parameter space where the SSPP has difficulty distinguishing dwarf and giant stars, and 
also constrain how many dwarfs may be lost with a stringent cut in logg. Finally, we use the Mg 



index to distinguish between giants and dwarfs in our sample ([Morrison et al.l 120031 ). 



For each line of sight in the SEGUE sample, we model the distribution of stars using TRI- 
LEGAL 1.4 (S I6.4.2p . We examine every target in these model distributions with — r) in the 
G/K dwarf range and log (7 <4.2 to determine the size of possible contaminatioqf|. Combining 
the proportion of subgiants along the TRILEGAL line of sight with the uncertainty in the SSPP 
surface-gravity measurements, we find that, over all lines of sight, the SEGUE color and magnitude 
criteria will include a mean of 2ibl% subgiants incorrectly identified as dwarfs. 

We then combine our population analysis with a study of the SSPP surface-gravity determi- 



*The surface-gravity limi t for these models is slightly higher than that applied to our sample, as these models are 
based on Padova isochrones l|Girardi et al.ll2004l ') , which have the main-sequence turnoff at log g k.A.2 whereas YREC 
isochrones place it at 4.1 
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nation using synthetic spectra. As the empirically-corrected YREC isochrones used in this study 
do not extend beyond the main-sequence turnoff, we adopt atmospheri c parameters for t he gi- 
ant and dwarf models from the Dartmouth Stellar Evolution isochrones (jPotter et al.ll2008l ) over 
a range of metallicity (see ? IB.4I for more details about the isochrones). Breaking the isochrones 
down into 0.01 magnitude blocks in color, we manufacture synthetic spectra fo r both giants and 



dwar f s using MARCS mod el atmospheres processed through TurboSpectrum (jGustafsson et al 



2008 



Alvarez fc Pleall998l ). O ur MARCS model atm ospheres assume solar-scaled abundances, 



with the solar composition from iGrevesse et al.l (|2007l ). and plane-parallel geometry. They cover 
a range in effective temperature from 4700-5800 K, —2.5 to +0.5 in [Fe/H], and 1.2 to 5 in log 5. 
Each of these spectra was adjusted to a range of S/N: 50, 25, and 10. We also analyzed the 
non-degr aded spectra. Th e spec tral synthesis and noise-modeling processes are discussed in more 
detail in lSchlesinger et al.l (j2010l ). 



Each of these simulated spectra is then processed through the SSPP using the parameter 
estimates as described in AppendixlAl A comparison of the calculated parameters and those input 
in the models for S/N '^25, the most common S/N for the G- and K-dwarf sample, is shown in 
Figure[6l At a S/N of 25, the surface gravity tends to be underestimated by 6%, less than the 
listed SSPP uncertainties. No giants or subgiants will be identified as dwarfs by their logg values. 
As S/N decreases, the uncertainties in the SSPP parameters increase. At the lowest signal to noise 
ratio in our sample, there is a 2% chance of a subgiant being identified as a dwarf by the SSPP 
and included in our G- and K-dwarf sample. From our population studies of the galaxy models, 
we expect approximately 2% of stars that fall into our G- and K-dwarf color cuts to be subgiants, 
and ~2% of these will be misidentified as dwarfs by the SSPP, resulting in a less than 1% chance 
that a subgiant will be counted as a dwarf in our sample. 



Our final check on the extent of evolved-star contamination uses the Mg index (jMorrison et al 



2OO3I ). At a given [Fe/H], giants will have a smaller Mg index value (i.e., less atomic Mg and 
MgH absorption) than dwarfs. As SEGUE targets K-giant stars in addition to K dwarfs, we focus 
on this spectral type, calculating the Mg index directly from the SEGUE spectrum. While this 
index has been calibrated for giants using known open and globular cluster members (Ma et al, in 
preparation), we have almost no calibrating observations of known metal-poor dwarfs. Therefore, 
we identified star s from our sample whic h were very likely to be dwarfs using a reduced proper- 
motion criterion (iLepine &: GaidosI I2OIII ). Binning in metallicity, we isolate a control sample of 
"true" dwarf stars by specifying that the total proper motion must be greater than 20 mas/yr, and 
the ro magnitude reduced proper motion is greater than 13. We then examine the spectra of targets 
where the Mg index is lower than expected for the true dwarfs, e.g., they are in the parameter space 
of evolved stars. This visual inspection consisted of three tests: 



The strength of the Mg feature in non-continuum-corrected SEGUE spectra (jMorrison et al 



2OOOI). 



A comparison of the strength of Ca I A4227 and Ca II K feature (jMorrison et al.ll2000l ). 
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The ratio of the strength of Sr II A4077 to three nearby Fe I hnes ([Morgan et al.lll943l : [Rose 



mi)- 



These three tests are valuable luminosity discriminants and isolate subgiant and giant stars mas- 
querading as dwarfs in the sample. Our visual analysis indicates that less than 1% of the total 
sample of SEGUE K dwarfs are actually evolved stars, confirming that subgiant and giant con- 
taminants have a negligible presence in our sample for both spectral types. We do not take this 
negligible contamination into account. 

Our analysis of synthetic spectra indicates that the SSPP tends to underestimate surface 
gravity, which prevents subgiants from entering the sample. However, it also means that a number of 
dwarf stars will fall out of the sample when it is selected using a surface-gravity cut. In particular, we 
will preferentially lose high-metallicity stars, as these have surface gravities closer to the boundary 
of logs'. We run a Monte Carlo analysis to examine how uncertainties in surface gravity will 
affect the MDF. We model each SEGUE hue of sight with the TRILEGAL Galaxy model fS [6X2]l . 
selecting G- and K-dwarf "spectroscopic" targets based on SEGUE photometric criteria. We then 
determine the "true" MDF for this sample. Next, we vary the modeled surface gravity values, 
convolving them with a Gaussian error distribution with a of 0.6 dex, the largest possible SSPP 
uncertainty. We then compare the MDF of these modified samples to the original, to estimate 
the MDF uncertainty from SSPP log 5 errors. For [Fe/H]> —1.0, the change in each MDF bin is 
around 3%. Below this metallicity, the percent change is higher, around 10%, due to the small- 
number statistics of this portion of the simulated distribution (Fig. [5]). These are much smaller 
than the expected MDF uncertainties from bootstrapping alone (^ 15. 2p and do not induce a bias in 
the metallicity distribution. 



4.4. Extinction 



Although the lines of sight in SEGUE typically probe above and below the plane of the Galaxy, 
undergoing small amounts of extinction, this small amount of reddening can affect SEGUE target 
selection, which is based o n color and magnitud e. We determine the extinction per plate using the 
r-band extinction from the Schlegel et al. ( 1998 ) maps for every object listed as a staj^. All but 11 
of the plates (1888, 2052, 2179, 2300, 2334, 2335, 2623, 2669, 2679, 2680, and 2805) have extinction 
less than 0.5 mags in r and Fi{g — r) of less than 0.2 mag. Reddening cuts remove around 6% of 
the sample. 

Although many of the lines of sight have minimal extinction, there is still a danger that 
reddening will bias the G- and K-dwarf sample. Closer stars wi ll have underg o ne les s extinction 
than those farther away; using the extinction corrections from I Schlegel et al.l (| 19981 ) may make 



^These reddening values were extracted from the DR7 table PhotoObjAlI. 
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close-in, metal-rich, stars appear artificially blue, leading to inaccurate distance estimates and 
bumping stars in and out of the G- and K-dwarf SEGUE sample. This should be a minimal effect, 
as our stars have a minimum \Z\ of 0.2 kpc, beyond the expected scale height of 0.14 kpc for dust 
in the disk (jMendez Sz van Altenalll998l ). 



(2011 
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sample of SEGUE main-sequence turnoff stars by matching the isochrone-predict ed {g — r) color, 
found using the SSPP parameters, to the observed. Unlike the lSchlegel et al.l (|l998l ) values, this red- 
dening calculation does not assume that each star lies behind the full amount of line- of-sight dust. 
For th eir four lines of sight with ^{g — r) less tha n 0.2 mag, sini il ar to our own, the ICheng et al 



(2011 



isochrone extinction is comparable to the ISchlegel et al.l (|l998l ) values over the full dis- 
tance range (Cheng , private communication). Thus, in areas of generally low extinction, using 
Schlegel et al.l ( 19981 ) extinction estimates will not induce a bias in our sample against metal-rich 



stars by preferentially removing nearby stars. 

As with photometric errors and undetected binarity, extinction can also artificially shift stars 
in and out of the G- and K-dwarf color and magnitude range. Specifically, if extinction is overes- 
timated, bluer stars will be scattered in and red stars out of the color range. To constrain these 
effects, we simulate the dis tribution of redden ing values at a range of extinctions, using the main- 
sequence turnoff stars from lCheng et al.l (j201ll ). on TRILEGAL Galaxy models along each SEGUE 



line of sight. Modeling the extinction with respect to distance along each SEGUE line of sight, we 
examine how the G- and K-dwarf sample changes. Although reddening does shift stars in and out 
of the G- and K-dwarf color range, it does not preferentially remove stars from a certain metallicity 
space. The uncertainty per metallicity bin in the MDF due to extinction is typically around 2% 
(see Fig.[5]). 



4.5. Volume Completeness 

Our G- and K-dwarf sample covers a wide range of metallicities. Metal-rich stars are brighter 
than metal-poor ones for a given [g — r) color. Thus, the volume coverage of the sample varies 
with respect to metallicity for the same magnitude range. At the bright end, metal-rich stars will 
saturate; at the faint end, metal-poor stars are not sufficiently luminous to be observed by SEGUE. 
To ensure that our SEGUE sample and MDF is not biased with respect to chemistry due to the 
metallicity-dependent volume coverage, we specify distance limits for our sample. 

We examine the YREC isochrones for a metallicity of [Fe/H] of -1-0.4 and —3.0, extracting the 
r-band magnitudes for targets with {g — r) of 0.48, 0.55, and 0.75. Using the distance modulus, we 
calculate the maximum and minimum distance for both spectral types over a range of metallicity. 
For K dwarfs, the distance range for SEGUE volume completeness is from 1.19 to 1.84 kpc; for G 
dwarfs, the distance range is from 1.59 kpc to 2.29 kpc. Note that these distances utilize adjusted 
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r-magnitude limits for the G and K dwarfs, as described in S I4.7.2I In this paper, we refer to the 
distance range from 1.59 to 1.84 kpc as the overlap distance range. These distance limits are used 
in comparisons between the G- and K-dwarf sample. For analysis of an individual spectral type, 
we use the larger distance range associated with either G or K dwarfs. We refer to this range as 
the spectral-type distance range. 

Random and systematic uncertainties in distance will shift stars in and out of the overlap range. 
We estimate the effect of these changes in distance on our MDF with a Monte Carlo analysis of the 
TRILEGAL simulation (^ I6.4.2p . Each star is assigned a typical percent change in distance from 
random and systematic effects, which is then convolved with a Gaussian over multiple iterations. 
We use a random distance uncertainty from errors in photometry, [a/Fe], and isochrone choice; the 
systematic uncertainties originate in age and binarity assumptions. For [Fe/H]> —2.0, variation 
in the systematic distance will change each MDF bin by around 3%. The number of stars at the 
metal-poor end is very small, leading to larger uncertainties, on average 35% (see Fig. [5]). The 
random distance uncertainties have a larger effect than systematic uncertainties, changing each bin 
by approximately 7% and 40% at the metal-rich and metal-poor ends respectively. These distance 
errors are included in our total MDF uncertainties. 

We do not include the change in distance due to SSPP [Fe/H] uncertainties as part of our 
random distance uncertainty, as these properties are strongly correlated. Overestimating the metal- 
licity will lead to an overestimate in distance, and vice versa. Uncertainties from [Fe/H] will both 
broaden the true metallicity distribution and shift stars in and out of the sample due to distance 
cuts. To estimate the size of these effects, we perform a Monte-Carlo analysis along the TRILEGAL- 
modeled lines of sight. For each modeled G- and K-dwarf star, we vary the [Fe/H], convolving the 
SSPP errors with a Gaussian, and adjust the corresponding distance to account for this change 
in metallicity. We then compare the different Monte Carlo iterations to the true distribution (see 
Fig- El)- The cumulative distributions of G and K dwarfs indicate that the [Fe/H] and distance error 
broadens the MDF of the TRILEGAL model. To quantify the extent of this effect, we estimate 
the slope of the cumulative distribution between fractions 0.25 and 0.75. Larger uncertainties in 
[Fe/H] will manifest in larger negative slopes in this portion of the cumulative distribution. For 
the true model distribution, the slope is around —0.6 dex~^, whereas it increases to —1.0 dex~^ 
when we factor in a^pe/n]- We list the broadened slopes for both spectral types in FigureEl Thus, 
uncertainties will cause us to measure a broader, by slightly less than a factor of 2, version of the 
true underlying distribution. 



4.6. Evolutionary Effects 

To derive a complete chemical history of the Galactic disk, we must ensure that no stars have 
evolved off of the main sequence, as this would bias our SEGUE sample against metal-r ich stars 



Figur e[2] shows the YREC isochrone set at an age of 10 Gyr for a range of metallicities (jAn et al 



20091 ). The G- and K-dwarf SEGUE color-cut regions are indicated. The most metal-rich 10 
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Gyr isochrones, calibrated to NGC 2862 and NGC 6791, have the main-sequence turnoff before 
reaching the red edge of the color range for the SEGUE G dwarfs. At younger ages, these metal- 
rich isochrones cover the entire G {g — r) range, indicating that there is a chance we are losing old, 
metal-rich G dwarfs due to evolutionary effects. 

To investigate the extent of this evolutionary effect, we compare the metal-rich samples of the 
two spectral types. Over the overlapping distance range for both G and K dwarfs, approximately 
3ibl% of G dwarfs have metallicities greater than 0. Around 10ib3% of K dwarfs are this metal rich. 
Between [Fe/H] of —0.3 and 0, the G- and K-dwarf distributions are consistent with one another. 
This indicates that our selection criteria loses metal-rich ([Fe/H]>0) G dwarfs older than 5 Gyr. 



4.7. Correcting for Metallicity Biases from Target Selection 

For our G- and K-dwarf sample to accurately reflect the chemical properties of these spectral 
types in the Galactic disk, we must ensure that we constrain and account for the different metallicity 
biases induced by SEGUE target selection. 

There are three metallicity biases in the SEGUE G and K dwarf sample that we correct for: 

• Other SEGUE target categories that are biased towards metal-poor stars overlap the color 
range of the SEGUE G and K dwarfs, biasing the color-selected sample in metallicity. 

• SEGUE has the same limited number of spectroscopic fibers for each line of sight, regardless 
of its stellar density. 

• The G- and K-dwarf color cuts select a range of stellar masses, and thus a subset of the mass 
function, that varies with metallicity. 

To correct for the effects of target selection, we compare the spectroscopic sample along each 
individual line of sight to the photometric sample. For each pointing and target category, we scale 
the spectroscopic distribution in tq such that it matches the underlying tq magnitude distribution 
of all of the stars of that could have been observed spectroscopically. Our sample consists of 152 
different lines of sight, each with a different distribution of spectral types and number of stars. 
Thus, we must correct for the target-selection biases on a plate-by-plate basis. 

For each bright and faint plate combination, we extract every possible photometric stellar tar- 
get within the plate radius that passes the G- and K-dwarf color and magnitude criteria, whether 
or not it was observed spectroscopicallfl We then match each of these stars with SEGUE spec- 
troscopy. Matching photometry and spectroscopy is a non-trivial task because some plates were 
observed multiple times and many stars were repeated, some unintentionally, as geometric overlaps 



®The photometric information is from PliotoObjAll database in DR7 
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from multiple plates, and others as purposeful re-observations. Additionally, targets from SEGUE 
were occasionally observed again as part of SEGUE-2. We removed all duplicate spectroscopic 
observations, selecting the highest S/N observations using the sciencePrimary parameter. We also 
trimmed our spectroscopic sample to eliminate any poor observations, removing all targets that 
have S/N <10, an incalculable [Fe/H], or specific warning flags (^ 12. 2p . These targets remain within 
the photometric sample, as they are simply counted as stars that do not have spectra. For each 
spectroscopic target that fulfills the G- or K-dwarf color and magnitude cut, we determine three 
weights: target-type, r-magnitude, and mass-function. As we do not know the log (7 of targets 
assigned SEGUE fibers in advance, we calculate these weights for all usable spectroscopic G and K 
stars. We trim our sample in surface gravity when we calculate our unbiased MDFs. 



4.7.1. Target-Type Weights 

The simple selection criteria of SEGUE G and K dwarfs overlap with the criteria for other 
different SEGUE targets. SEGUE categories often focus on specific ranges in parameter space, 
such as stars with low metallicity or small proper motions. Targets that fulfill multiple target- 
type criteria have multiple opportunities to be assigned a spectroscopic fiber. This leads to an 
overabundance of other stellar categories in the G- and K-dwarf sample. For example, the SEGUE 
low- metallicity star and K-giant categories use a ugr selection based around UV excess to identify 
targets, overlapping with the G and K dwarf criteria in {g — r). These target categories will bias 
the SEGUE G- and K-dwarf sample towards metal-poor stars (Figure[8|). 

Every stellar target observed photometrically by SDSS is labeled with the SEGUE categories 
it fulfills. Any individual star may pass the selection criteria for more than one target type and is 
then labeled as being a candidate for more than one SEGUE spectroscopic sample. We compare 
the number of stars identified as photometric candidates for each target type and each combination 
of target types (i.e., a K-dwarf and a low-metallicity target) with the number of those candidates 
that have good spectra, as defined in §[2j We assign a "target-type weight" such that the two 
distributions match. This removes biases due to the overlapping target types in SEGUE, be they 
over- or under-sampled by the spectroscopic fibers, ensuring that the types of stars probed by 
SEGUE reflects that of the underlying photometric sample. 

Figure[9] shows the mean target-type weight with respect to [Fe/H] for both spectral types. 
This correction reduces the proportion of metal-poor stars for both spectral types. Otherwise, 
the K-dwarf weights are typically much higher than those for G-dwarf stars. There is also more 
variation with [Fe/H]. While the selections of metal-poor and K-giant stars in the G-dwarf color 
region is relatively inefficient, the K-dwarf color selection is more contaminated with stars from 
other categories, biasing the metallicity structure. At [Fe/H]> —0.5, the K dwarfs have a target- 
type weight of around 2.2, indicating that metal-rich K dwarfs are under-represented by a factor of 
two. In contrast, the target-type weight for metal-rich G dwarfs is approximately 0.7, suggesting 
that they are over-sampled in SEGUE. Below [Fe/H] of —1.5, the target- type weights are around 
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0.7 and 0.4, for K and G dwarfs respectively. Our calculated target-type weights will be available 
with SDSS Data Release 9. 



Lee et al.l (j2011bl ) , I Bow et al.l (j2011al ) , and iLiu &: van de VenI (j2012l ) determine their stellar 



sample by extracting only stars assigned spectroscopic fibers as G-dwarfs. This greatly diminishes 
the size of the spectroscopic sample. Furthermore, although this method does not explicitly include 
stars targeted as low-metallicity, that does not make it unbiased in [Fe/H]. For example, if all of 
the low-metallicity stars along a particular line of sight receive SEGUE fibers as low-metallicity 
targets, none will be included in the G-dwarf sample, regardless of whether or not they fulfill the 
criteria, creating a bias towards metal-rich stars. Fortunately, the G-dwarf sample is significantly 
less affected by the metallicity bias of SEGUE target selection (Figure[9|). 



4-7.2. r-magnitude Weights 

SEGUE probes each line of sight with a limited number of fibers; the survey cannot spectro- 
scopically observe every target in a field. Lower latitude pointings are closer to the plane of the 
Galaxy, they have many more stars and tend to be more metal rich than high latitude lines of 
sight. SEGUE assigns the same number of fibers to each line of sight, regardless of latitude and 
stellar density. This leads to a bias against metal-rich stars in the observed MDF of the full survey. 
Furthermore, to accurately measure the properties of the Milky Way, we must ensure that our 
spectroscopic sample refiects the apparent-magnitude distribution of the underlying photometry 
along each line of sight. 

To calculate these "r-magnitude weights", we separate the spectroscopic and photometric 
sample into G and K dwarfs based on (g — r)o color. We examine the distribution in rp magnitude 
for each individual spectral type, binning up the spectroscopic and photometric targets in 0.5 
magnitude bins from tq =13 to 21.5 magnitudes. These bin sizes ensure that we typically have 
at least one spectroscopic target associated with each group of photometric stars. While the r- 
magnitude weights change when specifying a smaller bin size in tq, the metallicity distribution 
function changes very little, well within the expected uncertainties for each bin of [Fe/H]. We 
compare the number of spectroscopic targets in each bin to the number of potential targets in the 
photometric sample. The inverse of this ratio is the weight assigned to each spectroscopic target 
with a measurable [Fe/H] and S/N >10 to recreate the parent photometric distribution. This 
quantity, which we refer to as the r-magnitude weight, accounts for the fact that SEGUE does not 
have unlimited fibers. 

Examining the faint end of the G-dwarf sample reveals that there are very few usable spectro- 
scopic targets with S/N >10 at the faintest magnitudes. However, there are often many photometric 
targets at these faint magnitudes, leading to large r-magnitude weights and anomalous spikes in 
the metallicity distribution function. To determine the true magnitude range of the usable spec- 
troscopic sample, we examine the cumulative distribution of spectroscopic targets in rg for each 



-22- 



plate. We find that 85% of all spectroscopic G dwarfs have rg <18.45. By setting 18.45 as our 
faint magnitude limit for G dwarfs, we avoid anomalous weighting. The magnitude criteria for K 
dwarfs in SEGUE are significantly more conservative than that for the G dwarfs, extending to ro 
of 19; thus, the K-dwarf sample does not have the same faint-end problem. 

As K-dwarf stars are more frequent than G-dwarf stars, and sampled more sparsely in SEGUE, 
these targets tend to have higher r-magnitude weights. Both G- and K-dwarf stars also have high 
r-magnitude weights for lines of sights that probe lower Galactic latitudes. For G dwarfs with 
[Fe/H]> —0.5, the typical r-magnitude weight is around 13; the typical weight for K dwarfs is 
around 15. This value decreases to 8 and 11 for G and K dwarfs below [Fe/H] of —1.5, respectively. 
As with the target-type weights, the r-magnitude weights will be available to the public as part of 
SDSS Data Release 9. 



J^.1.3. Mass Function Weights 

The {g — r)o cut that defines the G- and K-dwarf targets corresponds to different mass ranges 
for each metallicity. For example, the G-dwarf color cut isolates stars with mass near O.7M0for 
[Fe/H] of —1.0 and near-solar mass stars for solar metallicity. Mass functions for the Galaxy predict 
a larger number of less-massive stars; t hus, the color cu t of SEGUE will result in a metallicity bias 



in our sample. Previous studies (e.g., I J0rgensenl |2000| ) discarded stars to find a mass range over 



which their samples were reasonably complete, but, due to our narrow color range, this would 
drastically limit our sample. The masses probed by our color and metallicity range are typically 
from 0.5 to 1.0 M0, but do not necessarily overlap for each metallicity bin. As this is not a large 
range in mass, we expect the metallicity bias introduced by our color cuts to be a small effect 
compared to other observational biases that affect the sample. 

We employ the TRILEGAL 1.4 model (^ I6.4.2p . which utilizes a Chabrier mass function, to 
estimate this effect. We extract two samples from the model for each individual line of sight. The 
first sample has the same color and magnitude range as the G- or K-dwarf sample and is restricted 
to stars with log (7 >4.1, to identify main-sequence stars. We bin this sample in metallicity for each 
spectral type and examine the distance range the sample covers at each [Fe/H]. The second sample 
is also limited to stars on the main sequence and in the G or K magnitude range. However, rather 
than a color cut, we extract all stars with masses between 0.5 to O.6M0that fall within the same 
distance range of stars in the color cut. This ensures we are comparing the two samples over the 
same volume of space. For each metallicity bin, we compare the number of stars that fulfill the color 
criteria to the number of targets with masses between 0.5 and 0.6 Mq within the same volume. This 
ratio indicates how much to scale each metallicity bin of the G or K dwarfs to simulate sampling 
a consistent portion of the mass function over the entire metallicity range. As the mass function 
varies slightly from plate to plate, we calculate these weights for each SEGUE pointing. 



We also calculate mass weights based on the ISchonrich Sz Binneyl (j2009al jbl) Galaxy models 
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(SB), which utihze a Salpeter IMF (jSalpeteijll955l ). to examine how much these values will vary 
with different assumptions about the mass function (S I6.4.ip . The distribution of weights from the 
SB models is comparable to the mass-function weights derived from the TRILEGAL simulation; 
using SB values rather than TRILEGAL has a negligible effect on t he metallicity distrib u tion. T he 
TRILEGAL models cover a wider metallicity range than those of Schonrich Sz Binnev ( 2009al fbl). 
reaching lower metallicities, so we use these weights throughout our analysis. 

As with the other two weights, the mass- function weights decrease the proportion of metal-poor 
stars. However, in contrast to the weights for overlapping categories and spectroscopic sampling, 
the mass- function weights are not large; above [Fe/H] of —0.5, both G and K dwarfs have a typical 
mass weight of around 1.0. For metal-poor stars ([Fe/II]< —1.5), the typical mass weight is 0.7 for 
G dwarfs and 0.9 for K dwarfs. These values will not be released as part of SDSS Data Release 9, 
as other groups will have their own preferred mass function. 



4-7.4- Testing the Target- Selection Weighting on Galaxy Models 

FigurefTOl shows the original and adjusted MDF for G and K dwarfs over all the SEGUE lines 
of sight for the spectral-type distance range. As explained in ^ l4.7.Hl¥.7.2l and l4.7.3l the corrections 
for SEGUE target-selection biases significantly affect the MDF, increasing the metal-rich end. 

To ensure that our weighting algorithm appropriately removes biases from target selection, 
we test it on a series of modeled lines of sight from the TRILEGAL 1.4 theoretical Galaxy model 
§ 16.4.21 We analyze each TRILEGAL line of sight as a photometric sample for SEGUE target 
selection, labeling all modeled stars which meet the color and magnitude criteria of the SEGUE G- 
and K-dwarf categories. We randomly select 375 G and 95 K dwarfs for each line of sight as our 
"spectroscopic" observations, weighting this mock SEGUE sample according to the methodology 
discussed in ^4.7.21 and ^4.7.31 We then compare the weighted parameter distribution in [Fe/H] 
in two ways to that of the underlying sample. First, we compare the "spectroscopic" sample, 
weighted in target type and r-magnitude, to the modeled photometric sample (see Fig. llip . The 
combination of these two weights recreates the metallicity distribution for all modeled stars that 
fulfill the G- and K-dwarf criteria. This suggests that our photometry-based weighting techniques 
allow us to examine the metallicity structure of the underlying stellar populations in the disk using 
our spectroscopic sample. 

Our second comparison, shown in Fig. 1121 uses the target-type, r-magnitude, and mass-function 
weight to recreate the distribution of all stars that meet the G- and K-dwarf magnitude criteria and 
fall in a uniform mass range. The underlying photometric sample that we tested the target-type and 
r-magnitude weights on is slightly biased in mass, as explained in § 14.7.31 The top plot compares a 
{g — r) color cut to a mass cut in metallicity space, where both samples cover the same volume for 
each individual metallicity bin. There is little difference between the two distributions. The bottom 
plot indicates that the color cut manifests as a bias towards metal-poor stars. Our target-selection 
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weights account for this, removing the color-cut effect such that it appears our sample is extracted 
from a uniform portion of the mass function. We also tested our target-selection corrections on the 
Galaxy model of Schonrich & Binney (SB, 2009a, b, S I6.4.ip . with similar results. 



4.8. Summary of the G- and K-dwarf Sample 

Our SEGUE sample consists of 16,847 K dwarfs with S/N >10, surface gravities larger 
than log 5 =4.1, appropriate ro magnitudes, and E((7 — r) less than 0.2 mag. This sample has 
-3.4<[Fe/H]<0.6, 4.1< log 5 <5.0, and 4150<Te// <6315 K. It contains stars with distances from 
0.3 to 5.8 kpc and covers Galactic radii from 4.8 to 13.4 kpc and heights from the plane up to ±4.0 
kpc. We have 24,270 G dwarfs with the appropriate S/N, surface gravity, extinction, and tq. This 
sample ranges from [Fe/H] of —3.1 to 0.6, surface gravities up to 5.0, and effective temperatures 
from 4210 to 6300 K. It extends to larger distances than the K-dwarf sample, from ~0.4 to 10.2 kpc, 
covering from 3.7 to 16.5 kpc in Galactic radius and -7.7 to 9.7 kpc in height above the Galactic 
plane. 

We have performed the first complete analysis of the different possible biases in this SEGUE 
sample. The color and magnitude selection criteria of the two categories result in a sample that 
is quite clean, with few stars of inappropriate color being shifted into the sample by photometric 
errors or undetected binarity. We also remove lines of sight from our sam ple with 'E{g — r) greater 



than 0.2 mag and establish that uncertainties in the lSchlegel et al.l (j 19981 ) extinction measurements 
have little effect on our distance estimates. 

We also constrain how stars of different evolutionary phase affect our sample. Combining an 
analysis of logs' with MgH and Galaxy models, we determine that the SEGUE G and K dwarfs will 
have a negligible amount of contamination from evolved stars. We will lose old, metal-rich stars 
from the G-dwarf sample due to evolution; our {g — r) color cut biases the G-dwarf sample against 
stars with [Fe/H]>0. 

We estimated the effect of different uncertainties in the G- and K-dwarf sample on the MDF 
(see Fig. [5]). The uncertainties are dominated by distance errors, both random and systematic. 
None of the various uncertainties, from photometric errors to undetected binarity, induce a bias in 
the metallicity of the sample. In our MDF, we include every different source of uncertainty in the 
structure, combining them in quadrature (^ 15. 2p . 

Although issues ranging from surface gravity to extinction have little effect on our derived 
MDF, the target selection algorithm of SEGUE does bias the sample in metallicity space, favoring 
metal-poor stars. The same number of stars are observed in fields at low Galactic latitude, which 
have a high stellar density and are dominated by high-metallicity stars, as in those at high Galactic 
latitude, where the stellar density is much lower and the stars are primarily metal poor. Thanks to 
the quantitative target-selection algorithm of SEGUE, we have both quantified and corrected for 
these different biases using a series of weights. The uncertainties from these values are integrated 
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into our MDF uncertainties via a bootstrap analysis (^ 15.2 



5. The Metallicity Distribution Functions of G and K Dwarfs over the Disk 

5.1. MDF of the Disk as a Whole 

For each spectral type, we determine the distribution of stars with respect to metallicity and 
then apply our target selection weights such that the distribution reflects the chemical structure of 
the underlying disk population. This allows us to improve on previous analyses, namely, photomet- 
ric metallicity estimates, limited volume coverage, observational biases, and/or small sample size. 
Figures[TU] and[T3] show the total MDF for G and K dwarfs over the spectral type and overlapping 
distance ranges, respectively. This is the first unbiased MDF using spectroscopic metallicities over 
a substantial volume of the disk for both G and K dwarfs. 

These MDFs provide information about the chemical structure of the disk from approximately 
0.2 to 2.3 kpc above and below the plane of the Galaxy, and from Galactic radii of around 6 to 
11 kpc (Figure fH]) . The MDF of both G and K dwarfs range from metallicities of around —2.5 to 
+0.4. The distributions peak at around [Fe/H] of —0.4 to —0.6, indicating that the sample has 
a large population o f thick-disk stars, in addition to substantial contributions froin the thin disk 
with [Fe/H]= -0.2 (|Gilmore Wvselll985l : Twvse fc Giln^ll995l : Ichiba fc Beer^ boool l . Unlike 
the dwarf stars in the GCS and RAVE surveys, the SEGUE sample has coverage of both the in-situ 
thin and thick disks. Thus, our total MDFs reflect the overall disk properties, rather than an 
individual component of the Milky Way disk. Detailed information about these distributions for 
further analysis is available in TableO 



5.2. Uncertainties in the MetalUcity Distributions 

We use a two-pronged approach to estimate the uncertainties in our MDF. First, we run a 
bootstrap analysis on our data set, which provides information about the variation in the sample, 
with respect to atmospheric parameters, photometry, and target-selection weights while taking 
into account that we do not know the true underlying metallicity distribution. We randomly select 
G and K dwarfs from the original spectroscopic sample until our modified version has the same 
number of stars of each spectral type, sampling with replacement, e.g., an individual target can 
show up multiple times in the bootstrap sample. We then examine the variation in the metallicity 
distributions of these iterations, comparing with our actual result. As the lower-metallicity bins 
have a smaller number of stars, their uncertainties are typically larger (Fig. [5]); below [Fe/H] of — 1, 
the uncertainty in the fraction in each bin is around 40%. Above this metallicity, the uncertainty 
in each bin is typically around 8%. 

Errors from photometry, undetected binarity, log 5, and extinction typically manifest as an 
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additional uncertainty per bin of around 5% each and induce no significant metallicity biases in 
the distribution (§[!]). However, random and systematic distance errors contribute a significant 
amount of uncertainty to our MDFs, around 40% for low-metalhcity bins (§[32]) ■ The total MDF 
uncertainties are a combination of the measured variation from all of these different sources, com- 
bined with those from our bootstrap analysis. These account for both the parameter and selection 
uncertainties and the fact that we do not know the true underlying distribution. We combine the 
different uncertainties in each MDF metallicity bin in quadrature; they are listed in TableO for 
both G and K dwarfs over different distance ranges. As the lower-metallicity bins have a smaller 
number of stars, their uncertainties are typically larger; below [Fe/H] of —2, the uncertainty in 
the fraction in each bin is around 70%. Between [Fe/H] of —2 and —1, the uncertainties drop to 
24%; above this metallicity, the uncertainty in each bin is typically around 10%. In addition, the 
correlated [Fe/H] and distance uncertainty will broaden the overall derived MDF (§1131 Fig. [7]). 



5.3. Variation with Height above the Plane 

The chemical structure of the disk in Galactic height is a useful diagnostic for dynamical 
formation models of the Milky Way. It also reveals the interplay of the proposed thin- and thick- 
disk by probing both components with an in-situ spectroscopic sample. 

Our G- and K-dwarf samples probe far above and below the Galactic plane; Figure[T5] shows 
the distribution of stars with respect to [Fe/H] and distance from the plane, using the spectral-type 
distance ranges for G and K dwarfs. We bin the sample in \Z\ and [Fe/H], adjusting each bin 
with the target-selection weights for each individual star. The colors represent the logarithm of 
the weighted number of stars in each bin. We compare the MDFs of the G and K dwarfs over the 
overlapping distance range at different heights above the plane of the Galaxy in FigurefTHl 



6. Discussion 



6.1. The MDF for Different Spectral Types 



Past work by iFavata et al.l (jl997l ) found a distribution of nearby K-dwarf stars that was signif- 
icantly more deficient in metal-poor stars than their G-dwarf sample. They conjectured that this 
result could originate from either a met allicity bias in their orig i nal ca talog or metal-enhanced cool 
star formation. However, later work bv iRocha-Pinto fc Maciell (119981) foun d their K-dwarf sample 



Kotoneva et al.l (120021 ') found a much broader 



(119971 ). Likely, the narrow dis- 



in good agreement with their G dwarfs. Similarly, 
distribution of K dwarf s with respect to [F e/H] than iFavata et al. 
tribution of K dwarfs in IFavata et al.l (j 19971 ) was due to small sample size and a tendency to target 
brighter stars for spectroscopy. 



We compare the MDF of the weighted SEGUE G and K dwarfs over the overlapping distance 
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range in Figure[T3l Our analysis confirms the work of iRocha-Pinto fc Maciell (|l998l ). which found 
the distribution of metalhcity in G and K dwarfs in the solar neighborhood to be in good agreement. 
For the overlap region for the G and K dwarfs, i.e., between 1.59 and 1.84 kpc, both spectral types 
peak at arou nd [Fe/H| of —0.5 , consistent with the peak metalhcity for the a-enhanced sample 



identified by iLee et al.l (|2011bl ). The peak agreement of the two spectral types implies that the 



chemical-enrichment and star-formation history for the G and K dwarfs are quite similar. 

The MDFs with respect to \Z\ of the two spectral types differ at the extremes of the metalhcity 
distribution (Figure [T6|) . For heights below 1 kpc, the K dwarfs have a larger population of metal- 
rich stars. As discussed in §1121 our color criteria for G dwarfs will remove the most metal-rich 
stars with ages greater than 5 Gyr. Above 1 kpc, there are few stars at or above solar metalhcity, 
so the discrepancy between the two spectral types is no longer noticeable. 

When we examine the metal-poor end for the G- and K-dwarf samples at different \Z\, we 
find a discrepancy (Figure [T7|) . For \Z\ distances below 1.5 kpc, the fractional contribution of stars 
with [Fe/H]< —1.0 agree within la (Figure[T7|) . The percentage of stars below [Fe/H] of —1.0 is 
consistent for the two spectral types up to \Z\ =1.5 kpc. However, at \Z\ greater than 1.5 kpc, 
the percentage of metal-poor ([Fe/H]< —1.0) K dwarfs is significantly larger than the metal-poor 
G dwarfs. For \Z\ >1.5 kpc, the fraction of G dwarfs increases from ~6% at smaller \Z\ to 9'^'^%. 
The K-dwarf percentage increases significantly more, from around 15% for 0.2< \Z\ <1.5 kpc to 
32l|2% above \Z\ of 1.5 kpc (see TablesHl [5l El [7|). For smaller \Z\, the median metalhcity for 
both G and K dwarfs agree to within 0.04 dex. For the highest \Z\ bin, the median metalhcity for 
the G dwarfs is —0.55, whereas the K-dwarf median [Fe/H] is —0.74, much more metal-poor. 

The origins of this discrepancy are currently unclear. The optimized version of the SSPP en- 
sures that the [Fe/H] estimates for the G and K dwarfs are consistent with one another; there is no 
systematic offset in SSPP metalhcity between the two spectral types of this size (Appendix lA.2|) . 
One possible cause of the MDF differences at large |^| is our assumed binarity fraction. If there is 
a larger probability of an undetected companion for K dwarfs than G dwarfs, we will be underesti- 
mating the distance to more K stars. This will cause more metal-poor stars to fall into the K-dwarf 
sample. 

Examination of the G- and K-dwarf MDF at different heights above the Galactic plane reveals 
that the K dwarfs are more metal-rich at low heights and more metal-poor at large \Z\. This makes 
the total K-dwarf MDF appear broader than the G-dwarf MDF over the overlapping distance range. 
As the fraction below [Fe/H] of —1.0 from both spectral types is consistent within 2a, we leave a 
more in-depth study of the metal-poor end of the G- and K-dwarf MDF for a future paper. 



6.2. Extending the G-Dwarf Problem 



The original G-dwarf problem stemmed from the observational results that a simple closed 
box model of chemical evolution over-predicted the number of metal-poor stars. In the top row of 
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Figure fTHt we compare our t otal MDF of G and K dwarfs to a simple closed box model with instan- 



taneous recycling, based on iPagell (|1997| ). As expected, we find that our empirical distributions 
of both G and K dwarfs are deficient at the low-metallicity end when compared to these models, 
confirming that the G-dwarf problem persists in the large Galactic volume of our survey and to 
lower temperatures. For the total sample, the fraction of K stars with [Fe/H] less than —0.8 is 
~17l;3%, whereas the model predicts 28%. Similarly, the model predicts 27% of G dwarfs should 



have [Fe/H]< —0.8, but we observe ~10 



+3c 



70. 



The G-dwarf problem was originally reported for stars within the solar neighborhood. Using 
the SEGUE sample, we find that the discrepancy with the simple closed box model is evident 
throughout the disk (Figure [T8|) . The discrepancies are on the order of 10-20% between the observed 
distribution and the simple closed box model, enhanced at larger distances from the plane, and 
consistent between the spectral types. 



6.3. Comparison of the G and K Dvi^arf MDF with Previous Work 



Our total MDF for G and K dwarfs covers a substantial portion of both the thin and thick 
in-situ Milky Way disk. As a conseq uence, the peak metallicity for both spect r al types is lower 



than previ ous analyses, such as RAVE (ISiebert et al. 



2011 



Co§kunoglu et al 



2011 



Bilir et al.ll2012l ) 



and GCS (jHolmberg et al.l 120071 : ICasagrande et al.ll201ll ). which are domin ated by thin-disk, more 
metal -rich stars. Similarly , the local MDFs for smaller samples, such as iRocha-Pinto &: Maciel 
(jl998l ) and iKotoneva et al.l (j2002l ). are more metal-rich than ours. These analyses have no stars in 
common with the SEGUE sample and probe a lower \Z\ range. Although the lowest \Z\ bin of G 
and K dwarfs ranges from 0.2 to 0.5 kpc, the median distance from the Galactic plane for these 
SEGUE stars is around 0.4 dex. Thus, our low \Z\ MDF is not directly comparable with local 
samples, and we defer this analysis to a later paper. 



Allende Prieto et al.l (j2006l ) measure the MDF of G-dwarf stars, identified by effective tem- 
perature, between \Z\ of 1 and 3 kpc. This volume of space is dominated by the thick disk. 



Allende Prieto et al.l ([200a) used the SDSS calibration star sample to examine the MDF between 
1 and 3 kpc, a volume of space dominated by the thick disk. The size and homogeneity of the data 
enabled measuremnts of the change in the observed sample MDF with respect to distance from 
the plane. However, this sample pre-dated the SEGUE survey, and the calibration star selection 
introduced significant metallicity and ev olutionary bias. The target selection of their sample causes 
the MDF of lAllende Prieto et al.l (|2006l ) to drop off steeply above [Fe/H] of —0.4; it may also cause 
the MDF to peak at around [Fe/H] of —0.7, more metal-poor than our MDF above \Z\ of 1 kpc. 



Whereas the GCS, RAVE, and lAllende Prieto et al.l (120061) sam ple probe individu al compo- 
nents of the disk, the photometric analyses of llvezic et al.l (|2008l ) and lBond et al.l (120101 ) examine 
a larger volume of space, from \Z\ of 0.8 to 7 kpc. At heights below 2 kpc, their MDF is slightly 
more metal-poor than ours, by around 0.1 to 0.2 dex, perhaps due to inaccuracies in their photo- 
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metric metallicity calibration. The contrasting MDFs of these different samples emphasizes that 
the chemical structure of the disk varies significantly with respect to spatial location. Only with 
a uniform, extensive, spectroscopic sample such as SEGUE can we obtain accurate constraints on 
the chemistry of the disk as a whole. 



6.3.1. \Katz et aU 1(201 i ) 



Past work by iKatz et al.l (|201ll ) sought to constrain the vertical metallicity gradient of the 
thick disk alone. Their sample consists of approximately 400 stars along two lines of sight; they 
calculate atmospheric parameters using the ETOILE spectral matching program. Utilizing a series 
of Besangon Galaxy models, they distinguish the thick-disk stars from the thin-disk population. 



The MDF of the iKatz et al.l (j201ll ) sample is slightly more metal-rich than our G- and K-dwarf 
sample at all \Z\ (Figure [T9|) . However, the tw o samples show ge neral similarities to one another 
in structure and peak metallicity. Although the iKatz et al.l (j201ll ) sample probes beyond the solar 



neighborhood, it is limited to only two lines of sight. Thus, it is promising for their analysis that 
we see general agreement in the metallicity distributions, but their limited sample size and volume 
make it difficult to measure the variation in the MDF at different locations in the disk. 



6.3.2. iLee et al\ 1(2011 a ) 



Lee et al.l (j2011bl ) use the [a/Fe] ratios of the SEGUE G-dwarf stars to chemically separate 
the thin- and thick-disk components. They use a different sample than ours, extracting only stars 
specifically targeted as G dwarfs with S/N >30, log 5 >4.2, and distances within 3 kpc. We 
compare our G-dwarf metallicity distribution over the spectral distance range to the a-poor and 
a-enhanced metallicity distributions in Figure[20l I n this figure, the two a-separated components 
are normalized by the total number of stars in the iLee et al.l (|2011bl ) sample. Farther away from 
the plane of the Galaxy, the a-rich component makes up a larger part of the sample, indicating 
that this a-rich population is an important contribution to our MDF far from the plane. 

At low h eights above the p lane, our observed MDFs contains a mix of the a-poor and a-rich 
samples fromlLee et al.l (j2011bl ). At \Z\ of 1.0 kpc, our MDF skews more metal-poor, matching the 
Lee et al.l (|2011bl ^ a-rich sample better. This suggests that, at these heights, the G- and K-dwarf 
sa mple is do r ninated by the a-enhanced disk stars. Based on this comparison with the results 
of iLee et al.l (l2011bl ) , our observed vertical variation in the MDF reflects the transition from a 
predominantly a-poor population at low |Z| to one dominated by a-rich stars at high \Z\, as well 
as a change in mean metallicity with distance from the plane. 
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6.4. Comparison with Galaxy Models 



Models of chemical and dynamical evolution have evolved substantially since the simple closed 
box model of ISchmidtl (|l963l ). Despite these improvements, models of Galaxy development remain 
hamstrung by the lack of extensive observational samples and are calibrated on samples from the 
local volume. Until SEGUE, there were few samples available to test and refine their predictions 
for the Milky Way properties far from the plane of the Galaxy. With our G- and K-dwarf MDFs, 
we can test existing Galaxy models and provide critical constraints on the metallicity structure 
with respect to location, to improve the model predictions and examine the star-formation history 
of the disk far beyond the local volume. In particular, the range of \Z\ heights makes the sample 
ideal for examining the transition between the thin and thick disk. 

We compare our MDFs with predictions from two Galaxy models. Those of Schonrich &: 
Binney (SB, 2009a, b) combine both chemical and dynamical evolution, whereas the TRILEGAL 
1.4 model is based on empirical relationships. The models differ importantly in their modeling of 
the thick disk. While the TRILEGAL model treats the thick disk as a separate component of the 
Galaxy and allows different input parameters for its spatial distribution and local density, the SB 
models implement the radial migration of disk stars, forming a thick disk in the solar neighborhood 
via outward migration of older, inner disk stars. Thus, the SB model's metallicity distribution at 
high \Z\ will be constrained by the chemistry of old, inner disk stars, while TRILEGAL 's thick 
disk metallicity distributi on is fixed. We a re also currently working on a comparison with the 
Besancon Galaxy model ([Robin et al.l 120031 ) (Schlesinger et al., in preparation). We also have 
provided nun ierical tables of ou r MDFs (Tables O [U O [U [7]) for comparison with future Galaxy 
models (e.g.. ICalura et al.ll2012l ). 



6.4- 1- The Schonrich & Binney Galaxy Models 



The models of Schonrich fc Binney (j2009al lbl) rely on assumed star-formation laws, nucle- 



osynthetic yields, and a Salpeter initial m ass function (ISalpetei 



19551). These models base th eir 



star-formation rate on the Kennicutt law ( Kennicutt et al. 1998 ): Schonrich &: Binnev ( 2009al lbl) 
model the expected nucleosynthetic yields in conjunction with inflows, resulting in a predicted 
age-metallicity relationship, rather than building off an empirical version. In addition, these 
models inclu de dynamical behavior, spec ifically, radial mixing, which alters the modeled chemi- 
ca l structure. ISchonrich &: Binnevi (l2009bl') compare the modeled vertical density structure to that 



of Juric et al. (2003) and 



Ivezic et al.l (120081 ) , finding that the model accurately predicts the density 



structure of various Gala xy components. Finally, these models are built using the BaSTI isochrones 
(jPietrinferni et al.ll2004l ). which they artificially widen to simulate the photometric errors of SDSS. 
Thus, the expected uncertainties derived from applying our techniques on these models should 
accurately reflect the photometric errors of the SEGUE sample. We do not factor in the other 
uncertainties in the SEGUE sample, such as [Fe/H] and distance, into the SB models. 
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The SB model h as been cahbrated in [Fe/H] with th e Geneva-Copenhagen sur v ey (INordstrom et al 



2004 



Havwoodl[2008l ^ and with updated parameters from lCasagrande et al.l (j201ll ). ISchonrich &: BinnevI 



(j2009al ) have shown that their model accurately reflects the properties of stars in the solar ne i ghbor - 
hood. They compar e their model predict ions t o other ob s erved samples, such as iFuhrmannI (|l998l ). 



Bensbv et al.l (j2003l ) . iReddv et al.l (j2006l ). and iHavwoodl (|2008l ). to constrain chemical parameters. 



Although many of these samples contain thick-disk stars, they are all limited to t he solar neighbor- 



hood , which is thin-disk dominated. The GCS sam ple is typically withi n 200 pc (jCasagrande et al 



20 111 ). Similarly, the farthest distance probed by the lReddv et al.l (j2006l ) sample, which was selected 



based on Hipparcos data, is 150 pc. Thus, although the SB Galaxy model predicts the structure 
of the disk as a whole, it is only calibrated on the local volume. 



6.4.2. The TRILEGAL Galaxy Models 



Whereas the SB models simulate the disk's chemical and dynamical evolution, the TRILEGAL 
1.4 models are more empirical in scope, utilizing a derived luminosity function, built from an 
assumed star-formation rate, age-metallicity relation, and initial mass function. These parameters 
vary slightly for different components of the Galaxy. 

The TRILEGAL 1.4 program allows us to simulate different lines of sight and assume a range 
of Galaxy parameters. We manufacture a TRILEGAL model for each SEGUE line of sight, using a 
limiting magnitude of 23.0 in r, no dust extinction, and no binary contamination over a circular 7 
square degree area of the sky. We specify an exponential thin disk with a two-step star-formation 
rate. The thin disk is assi gned a scale length of 2800 pc and a scale height which increases with 
stellar age, as described in lGirardi et al.l ( 20051 ) . The assu med age-metallicity relationship for this 
component is based on the empirical determination from iRocha-Pinto et al.l (|2000l ) , and includes 
a-enhancement at lower [M/H]. The thick disk is also modeled as an exponential, with a scale height 
of 800 pc and a scale length of 2800 pc. In contrast to the thin disk, a constant star-formation 
rate and a Z of 0.008 ([Fe/H]= —0.7 for an a-rich population) and cr[M/II] of 0.1 dex (which is 
typically smaller than metallicity measurement errors) is assumed for the thick disk. In addition to 
the thin and thick disk, we specify an oblate spheroid halo and no central bulge. These properties 
are combined with the Padova evolutionary tracks to produce a derived luminosity function. Note 
that the TRILEGAL models predicts the total metallicity, [M/H], rather than [Fe/H]. 

This model shows good agreement between the observed and predicted star count s when com- 
pared to various stellar catalogs, such as 2MASS and Hipparcos (jGirardi et al.ll2005l ). There has 
been little chemical calibration of the TRILEGAL model. 
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6.4-3. Accuracy of Galaxy Models 

We compare the G-dwarf metallicity distribution to the SB and TRILEGAL Galaxy models 
in Figures[2T] and [22l respectively. For both models, we combine all of the SEGUE-modeled lines of 
sight and extract all simulated stars that fulfill the color, magnitude, and distance criteria for our 
G-dwarf spectral type. 

Both of these models are calibrated to past analyses, which are typically limited to the local 
volume. Thus, we expect agreement between the model and SEGUE sample close to the plane of 
the Galaxy. However, Figures[2T] and[22] indicate that, although there is substantial overlap in the 
range of [Fe/H], both models are more metal-rich and sharply peaked than our observed sample. 
The observed sample is broadened significantly by correlated errors in [Fe/H] and distance, in 
addition to general uncertainites in [Fe/H] (see ^ 15.21 and Figure[7]). The models do not simulate 
these uncertainties and thus exhibit a narrower MDF. 

At higher distances from the plane, the metallicity discrepancy between the SB model and 
the SEGUE sample increases. Whereas the SEGUE distributions become more metal-poor with 
increasing |Z|, the models remain at an [Fe/H] of —0.1 for SB. In fact, the SB distributions become 
slightly more metal rich at greater distances from the plane. We find similar discrepancies for the 
K-d warf sample. We co mpare the SB modeled lines of sight with the a-separated stellar populations 



from lLee et al. ( 2011b ) in FigureEH The model agrees well with the a-poor distribution alone but 



are unable to accurately represent the a-rich component. 

Rather than [Fe/H], the TRILEGAL model predicts [M/H], which includes contributions from 
[a/Fe]. The modeled [M/H] will thus be systematically higher than our [Fe/H] at large \Z\, where 
[a/Fe] is expected to be large. We model t he relation s hip be tween [a/Fe] and [Fe/H] using a 



simple linear relationship based on the work of lLee et al.l (j2011bl ). and remove the expected [a/Fe 



contribution from the [M/H] value of each simulated TRILEGAL star. We assume that stars of 
solar metallicity have [a/Fe] of 0. As [Fe/H] decreases, [a/Fe] increases to 0.20 at [Fe/H]= —0.5, 
and 0.40 at [Fe/H]= —1.0. Without this correction, the distribution is similar to that of SB, 
where the model is systematically more metal rich than the observations at all heights. When 
we account for [a/Fe], the match between our observed G-dwarf MDF and the model improves, 
however, TRILEGAL is stiU offset in [Fe/H], by around 0.1-0.2 dex, at all heights (Figure[22]) . We 
find that the model does not adequately simulate the structure of the a-rich population at large 
heights. In particular , from 1.5 to 2.3 kpc in \Z\, the TRILEGAL distribution is very similar to 



that of the iLee et al.l (j2011bl ) a-poor component. Note that we can improve the match between 



the simulation and observations by adjusting the contribution from [a/Fe]. 

The two models predict stars at large \Z\ to be more metal rich than we observe. The current 
inability of these Galaxy models, in particular the SB model, to predict the chemical structure of the 
disk points to the danger of assuming that the local volume reflects the properties over the span of 
the disk. Past analyses of the metallicity distribution, used to calibrate these Galaxy models, were 
limited to the solar neighborhood; they were dominated by thin-disk stars. While these samples 
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contain thick-disk stars passing through the local volume, they may not provide a complete picture 
of the component's structure. There appears to be a significant disagreement between the picture 
of the disk structure from local kinematic samples and our in-situ observations. Thus, our unbiased 
metallicity distributions of the SEGUE G and K dwarfs, and how they change with respect to R 
and \Z\, provide critical constraints for refining and improving Galaxy models, such that we gain 
a better understanding of the different chemical and dynamical processes occuring in the disk and 
no longer rely on local samples alone. 



7. Summary 



In this work we examine the variation in [Fe/H] of the SEGUE G- and K-dwarf sample with 
respect to position in the Galaxy. These distributions allow us to better constrain the chemical 
and dynamical evolution processes in the disk. To ensure the accuracy of this sample, we optimize 
the SSPP to better handle low-temperature dwarfs and develop and test a procedure to account 
for target-selection biases in the SEGUE sample, using SDSS photometry for each individual line 
of sight. This technique corrects for the metal-poor bias in the spectroscopic sample and is critical 
for accurately constraining the metallicity structure of the Milky Way disk. For the first time, we 
can utilize an unbiased sample of both G and K dwarfs with accurate atmospheric parameters to 
examine the chemical structure over a large volume of the disk. 

Using distances estimated with an isochrone-matching technique, accurate to ~8% for metal- 
poor and ~18% for metal-rich stars, we determine the MDFs of the G- and K-dwarf samples 
(Fig. lisp and examine how the metallicity structure of the disk varies with respect to position and 
spectral type. Both the total MDFs and MDFs with respect to height above the plane of the Galaxy 
for the G- and K-dwarf samples are in good agreement with one another, implying that the two 
spectral types undergo similar chemodynamical evolution over the span of the disk. We compare 
both distributions to the simple closed box model, finding that the G-dwarf problem persists to 
both lower temperatures and is prese nt throughout t he volume of the Milky Way disk (Fig. 1181) . 
This expands on previous wo rk, su ch aslMouldl (119821) .iFavata et al.l (|l997l ). lFlvnn &: Morelll (|l997l ). 



Rocha-Pinto &: Maciell (|l998l ). and lKotoneva et al.l (|2002l ). which examined the solar neighborhood. 



This is the first work where we can confirm that this discrepancy extends past the solar neighbor- 
hood with a large and unbiased stellar sample of two cool spectral types. In particular, all previous 
analyses of K dwarfs were much more limited in sample size and volume. 

Analyses of the local volume find that the thin and thick disk are distinct with respect to 
[a/Fe]. As it is less likely than kinematics to change significantly over time, stellar chemistry 
is a useful diagnostic to disentangle the two components. Our analysis of the G- and K-dwarf 
MDF finds that both spectral types have more metal-poor stars as \Z\ increases (Fig. 1161 Tables^ 
ying similar star- formation histories. Comparison with the [a/Fe] analysis of 
suggests that our vertical metallicity structure represents the transition between 
a population composed largely of a-poor stars at low \Z\ to predominantly a-rich stars at high 
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\Z\. This change with respect to \Z\ is not predicted by the SB or TRILEGAL Galaxy models. 
Both are more metal-rich than our observed distribution. In particular, the SB model appears to 
simulate only the a-poor component of the Galaxy, suggesting that there are discrepancies with the 
picture of the extended disk presented by local samples and our own in-situ observations. Thus, 
our observed MDFs provide a valuable constraint on future Galaxy formation and evolution models 
and strongly warn of the dangers of using existing models to predict behavior far from the local 
volume. 
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A. Optimizing the SSPP for Cool Dwarf Stars 

The SSPP analyzes a wide range of spectral types, and not all of the techniques employed are 
appropriate for cool dwarfs. To ensure that the SSPP provides accurate atmospheric parameters 
for G and K dwarfs, we have developed an optimized version of the SSPP, carefully selecting, and 
occasionally adjusting, the log 5 and [Fe/H] techniques used in the analysis. 
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A.l. Surface Gravity Determinations 

The SSPP uses ten m ethods to estima te surface gravities for stars, each appropriate for a 



different range in {g — r)o (jLee et al.ll2008al ). Although all of the techniques cover the range of 
(g — r)o for G and K dwarfs, they are not all similarly accurate. As shown in Fig. l23l the CaI2 
and MgH indicators have large tails that extend to low values of log (7 not seen for other indicators. 
Comparison between log g estimates from Kepler and those from the SSPP indicate that the CaI2 
and MgH techniques return systematically low surface gravities. To confirm this, we tested the 
different surface gravity indicators on synthetic spectra (see § 14.31 for more information); while the 
majority of the methods accurately reported the input log g, both CaI2 and MgH exhibited false 
low surface gravity tails. Examination of the SSPP revealed that the implementation of these two 
techniques was flawed; both CaI2 and MgH utilized an inaccurate conversion between (g — r) and 
(B-V). Furthermore, we found that the relationship implen iented in the SSPP for these techniques 



was an outdated version, different than that reported in iMorrison et al.l (|2003l ) . In addition to 
removing CaI2 and MgH, we also eliminate WBG technique from our log 5 estimate. As it was 
designed for evolved stars, this indicator rarely works well for the G- and K-dwarf sample. 

Using the CaI2, MgH, and WBG to determine logs' will shift the calculated values lower, 
potentially eliminating some dwarfs from our sample by making their surface gravities indicative 
of evolved stars. For our targets, we use log 5 values which do not take these three methods 
into account. This tightens the distribution of surface gravities, diminishing the population that 
artificially lies in the low- valued tail (see Fig. l23p . On average, removing CaI2, MgH, and WBG 
increases the surface gravities from DR8 by less than 0.1 dex (see Figure[24 



A. 2. [Fe/H] Determinations 

As with the surface-gravity estimates, not all of the metallicity-estimation techniques available 
in the SSPP are accurate for G and K dwarfs. For example, the CanK2 and WBG distributions 
are erratic and spiky (Figure [25]) . Other techniques, such as CanK3, ACF, ANNSR, and CaHT, 
are limited in the metallicity range they cover, only reaching as high as solar metallicity. This limit 
will bias their measurements towards lower metallicities. We exclude these six techniques from our 
[Fe/H] determinations. 

Finally, there are other methods in the SSPP that appear to have well-behaved distributions in 
[Fe/H] but are inappropriate for our sample. Two of the Kurucz-model-based methods (KI13, K24) 
are designed for stars with temperatures greater than 5000 K; there are numerous K dwarfs for 
which this method is not accurate. The ANNRR method is also problematic; thi s neural-network 



techn ique is based on outdated parameters determined by the SSPP for DR7 (jAbazajian et al 



2009 1. Thus, our [Fe/H] estimates are based solely on NGS2, NGSl, and CaHKl. 



Although these three remaining methods cover the appropriate range of parameter space, they 
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still require adjustment. In particular, all three show artificial peaks. These methods utilize a 
minimization technique in Tefj, log 5, and [Fe/H] to match target spectra with various synthetic 
grids. These grids cover a fairly narrow range of parameter space, resulting in the sharply peaked 
structure. By mo difying the inte rpolation method in parameter space from the [Fe/H] values 



reported in DR8 (jAihara et al.ll201ll ). these peaks are smoothed, providing a more realistic distri- 
bution in [Fe/H]. Figure[24l indicates that, on average, these modifications have very little effect on 
the estimated [Fe/H]. These changes will be detailed in the upcoming SEGUE DR9 paper, and are 
outside the scope of this work. 



B. Uncertainties in Calculated Distances 
B.l. Photometry 

The photometric uncertainty of SDSS targets is typically 2-3% for each filter. Using the 
TRILEGAL simulations along SEGUE lines of sight, we simulate the photometric uncertainties in 
SDSS based on a Gaussian distribution. These changes factor in to distance uncertainties twice; 
first, when we match each target to an isochrone in {g — r)o and second, when we calculate the 
distance in each individual filter. The change in distance with respect to photometric uncertainty 
is typically less than 6% over the entire metallicity range. 



B.2. [Fe/H] 

There are two main components of uncertainty in the SSPP [Fe/H] estimates. First, there is 
the error which refiects the variation from the different techniques used by the SSPP to estimate 
[Fe/H]. This uncertainty from the three methods we utilize has a mean value of 0.07±0.06 dex. 
Second is the uncertainty determined by comparing with stars of known metallicity, from clusters 
and high-resolution observations. This uncer tainty varies with respect to S/N, ranging from 0.24 



dex for S/N of 50 to 0.45 dex for S/N of 10 (jLee et al.ll2n08al l 



For each SEGUE line of sight, we examine how the S/N varies with respect to magnitude, 
simulating the relationship for the associated TRILEGAL model. Based on the assigned S/N, 
we estimate the expected [Fe/H] uncertainty, combining it in quadrature with a simulation of the 
pipeline [Fe/H] variation, convolved with a Gaussian. We then compare the original distances calcu- 
lated for the TRILEGAL stars to those using [Fe/H] values adjusted for the simulated uncertainty. 
Changes in [Fe/H] have a larger affect on distance estimates for metal-rich stars ([Fe/H]> —1.0); 
these dominate the uncertainty with a distance error of 15%, whereas more metal-poor stars have 
an error of around 4%. 
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B.3. 



[a/Fe] 



An et al.l (j2009l ) assign each YREC isochrone a specific [a/Fe] value based upon its [Fe/H]. 



These values may not align with the [a/Fe] of the individual target, which will change the distance 
estimates, espe cially for the meta l-rich regime where small changes in chemistry significantly change 
the isochrone. iLee et alJ (|2011al ) found that for a S/N >20, the SSPP can determine [a/Fe] with 
an uncertainty of approximately 0.1 dex. At lower S/N, the uncertainty increases to 0.2 dex. 
To estimate the change in estimated distance associated with a 0.2 dex uncertainty in [a/Fe], we 
compare the magnitudes of the Dartmouth isochrones at each value of [Fe/H] with [a/Fe] of to 
those with [a/Fe] of +0.2 and —0.2. We must use the Dartmouth isochrones to estimate this effect, 
as the YREC isochrones do not model a range of [a/Fe] for each [Fe/H]. We simulate the change 
in Dartmouth ugriz at each [Fe/H] from a change in [a/Fe]; similar to the distance uncertainties 
from [Fe/H], the change in distance is larger, approximately 10%, at the metal-rich end. Below 
[Fe/H] of —1.0, the o'rf([a/_Fe]) is approximately zero. 



B.4. Isochrone Choice 



Each individual isochrone set will yield slightly different ugriz values at the same {g — r), 
resulting in systematic differences in distance with respect to one another. The various assumptions 
integrated into each of these different isochrone sets mea ns that for each {q — r) color, the re is a 
range of appropriate r magnitude about the YREC value (jCallart. Zoccali. &: Aparicidl2005l ). The 
difference between isochrone sets is more extreme at the metal-rich end, making the uncertainty 
due to isochrone choice larger than for metal-poor stars. 

Distances calculated using YREC isochrones tend to be smaller than those calculated via 
Dartmouth isochrones by around 6 % for stars with [F e/H]> —0.5; below this, the offset is around 



3%. In contrast, both the BaSTI (jCassisi et al.l 120061 ) and Padova (iGirardi et al.ll2004l ) isochrone 
sets will measure distances smaller than those from the YREC isochrones, approximately 8% for 
[Fe/H]>0. For more metal-poor stars, the differences between BaSTI and YREC are around 1%. 

Changing which isochrone set we use will induce a systematic change in our estimated dis- 
tances. However, when we consider the many different isochrones, choosing a particular set can 
make our distances smaller or larger. Thus, we treat isochrone choice as a random uncertainty of 
approximately 7% for high-metallicity stars ([Fe/H]> —0.5) and 3% for low-metallicity stars. 



B.5. Age 

To calculate distances, we assume all stars have an age of 10 Gyr and compare them to 
isochrones of the same age. However, the actual ages of these stars are unknown. In particular, 
the most metal-rich stars are likely to be significantly younger than this value, causing us to 
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systematically overestimate the distances to these stars. 

To estimate the uncertainty from age, we utilize the modeled lines of sight from the SB Galaxy 
models. These models simulate a wide range of ages, from 2.5 to 11.5 Gyr using an age-metallicity 
relationship built off of expected nucleosynthetic yields (SI6.4.1I). In contrast, th e TRILEGAL 
Galaxy model assumes an age-metallicity relationship from lRocha-Pinto et al.l ([20001), and contains 
stars from only 7 to 9 Gyr in age. For each modeled line of sight, we calculate the distances for 
simulated G and K dwarfs using 10 Gyr YREC isochrones. We then calculate the distances for 
the same sample using YREC isochrones that match the target in age, examining the differences 
between the two values. For the youngest stars, from 2.5 to 5 Gyr, using 10 Gyr isochrones 
creates a systematic 3% overestimate of distances. This value decreases as the ages of the stars 
increase. With respect to metallicity, this means that the metal-rich stars, with [Fe/H]> —0.4, will 
be systematically offset in distance from 1 to 3%. More metal-poor stars are unaffected. As all of 
our targets are on the main sequence, it is not surprising that we see little effect with respect to 
age assumptions; at this point in stellar evolution, isochrones of a given metallicity overlap quite 
closely. 



B.6. Binarity 



An undetected companion will make an individual target appear brighter and redder, causing 
us to underestimate the distance to the p r imary star. We simulate a population of binaries, using 
techniques described in ISchlesinger et al.l (j2010l ) , to examine how much undetected bin arity will 
affect our distance calculations. We model a population of binaries, using a Chabrier IMF (jChabrier 
20031 ) for both the primary and secondary stars, and determine the change in each ugriz filter with 
an undetected companion over a range of metallicity. We then match up each star in the TRILEGAL 
modeled lines of sight with these ugriz changes, comparing their original YREC distances to those 
estimated with the companion present. Binarity will typically lead to a systematic distance offset 
of 5%. Metal-poor stars ([Fe/II]< —1.0) have a slightly larger change in distance, around 7-8%, 
whereas more metal-rich stars are around 4%. Although small changes in magnitude have a larger 
effect on distance estimates for metal-rich stars, the changes in magnitude from binarity are larger 
for metal-poor stars, such that t he change in distance is slig htly larger. As approximately 65% of 
all G and K stars are in binaries (jPuquennoy fc Mayoiill99ll ). this results in a systematic offset to 
smaller distances of around 3%. 



In addition to affecting photometry, undetected pairs can also alter our distance measurements 
via changes in the metallicity deter mination. The expected [Fe/H] uncertainty at a S/N of 25 is 
±0.30 dex. ISchlesinger et al.l (|2010l ) find that around 95% of all pairs will have a shift in [Fe/H] 
less than 0.15 dex. Thus, distance uncertainties stemming from a secondary's effect on [Fe/H] 
determinations are well within the error expected from the SSPP [Fe/H] determination, and we do 
not take them into account separately. 
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B.7. Summary of Random and Systematic Uncertainties in Distance 

We have explored the different uncertainties that factor into our distance calculations using 
simulated lines of sight from TRILEGAL and SB Galaxy models. The random uncertainty in 
distance stems from photometry, [Fe/H] determinations, and [a/Fe] estimates. Photometry con- 
tributes a 6% uncertainty over the entire metallicity range; the chemistry of a star however, causes 
a different amount of uncertainty for metal-rich and metal-poor stars. Metal-rich stars will have a 
change in distance of 16% from [Fe/H] and 10% from [a/Fe] uncertainty. At the metal-poor end, 
these values are smaller, around 4% from [Fe/H] and 0% from [a/Fe] uncertainty. In addition, we 
could use a number of different isochrones to estimate our distances, adding an additional random 
uncertainty in distance of 10% for metal-rich and 3% for metal-poor stars. Combining these in 
quadrature, we estimate our random distance uncertainty ranging from 15—18% for [Fe/H]> —1.0 
to 8% for stars with [Fe/H] below -2.0 (Fig. ED . 

Our distance estimates also suffer from systematic uncertainties from age assumptions and 
binary contamination. While the distance uncertainties from age assumptions vary with metallicity, 
from to 3% from the metal-poor to metal-rich end, binarity results in a more consistent offset of 
around 5% over the full [Fe/H] range. 



B.8. Comparison with llvezic et al.l (|2008l ) Photometric Parallax 



Ivezic et al.l (|2008l ) develop a photometric paral lax relationship which takes into account metal- 
licity, building upon the work bv ljuric et al. (2008). The shape of the photometric relationship is 
based on SDSS observations of five globular clusters: M2, M3, M5, M13, and M 15. These clus- 
ters range from 7 to 12 kpc in d istance and —1.2 to —2.1 in [Fe/H] (iHarrisI 119961 ). Using clusters 
from IVandenberg &: ClemI ( 20031 ). which range from —2.5 to -|-0.12 in [Fe/H], they then constrain 
the extent of magnitude shifts for different metallicities. Combining the equations determined for 
different ranges in {g — i), their final relationship is as follows: 



Mr{g - i, [Fe/H]) = M^ig - i) + AM,([Fe/i/]) 



(Bl) 



M°(ff-i) = -5.06-^14.32(5- 



12.97{g-i)'^ +6.127 {g - 



1. 267 {g - r)'^ + .0967 (g-i)^ (B2) 



Mri[Fe/H]) = 4.50 - l.ll[Fe/H] - 0.18[Fe/if]^ 



(B3) 



Ivezic et al.l ()2008l ) specify that these equations are valid for dwarfs with 0.2<{g — z)<4.0. Less 
than 1% of G and K dwarfs fall outside of this color range. The y expect uncertai nties of around 
5% and 10-15% from systematic and random errors, respectively (jBond et al.ll2010l ). 



We compare our YREC isochrone distances to those calculated using photometric parallax, 
finding an offset of approximately 8%, with YREC distances typically smaller. This result is due 
to an ~0.2 mag offset between the YREC isochrones and photometric metallicity relations. This 
offset varies with respect to metallicity: for [Fe/H]>0, the change in distance is around 6%, whereas 
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for [Fe/H] between —1.5 and —2.0, the change is minimal. Finally, for the most metal-poor stars 
([Fe/H]< —2.0), the photometric parallax distances are actually smaller than those determined 
using YREC isochrones, by around 7%. To ensure our distance method does not significantly 
change our results, we have also calculated the metallicity distribution functions using photometric 
parallax distances instead of those from YREC. There are only slight changes, all of which are well 
within the expected bootstrap uncertainties for the distributions and gradients. 
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Fig. 1. — The atmospheric parameters of our sample of around 26,600 G and 18,500 K dwarfs with 
appropriate S/N and surface gravity. The sohd hne represents the G dwarfs, and the dotted hne 
the K dwarfs. All of these stars have logg >4.1, S/N >10 per pixel, where each pixel is ~lA, and 
a measured temperature and metallicity. 
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Fig. 2. — The YREC isochrone set for a range of metallicities at an age of 10 Gyr. Each isochrone 
is labeled with its specific [Fe/H]. The vertical dashed lines show the (g — r) color range for SEGUE 
G and K dwarfs. These isochrones extend slightly past the main-sequence turnoff to the subgiant 
branch. 
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Fig. 3. — The top plot shows the distribution of G- (sohd) and K-dwarf (dashed) distances, cal- 
culated by matching to YREC Isochrones. The K dwarfs tend to be closer in distance than the 
G dwarfs, which we expect as they are typically less luminous. The K dwarfs are peaked at a 
distance of around 0.95 kpc; the G dwarfs peak around 1.25 kpc. The overlapping distance range 
for K dwarfs and G dwarfs ranges from 1.19 to 1.84 kpc and 1.59 to 2.29 kpc, respectively. The 
middle plot shows the distance error due to random errors with respect to [Fe/H]. This includes 
uncertainty from SSPP [Fe/H], photometry, [a/Fe], and isochrone choice. We divide the sample 
into 25 bins of [Fe/H]; the points represent the mean uncertainty in each bin. The bottom figure 
shows the systematic uncertainty with respect to [Fe/H], plotted in the same way. This combines 
the systematic underestimate in distance from assuming all stars are alone and the overestimate 
from assigning an age of 10 Gyr to each star. Although we use TRILEGAL for the majority of 
these uncertainty determinations, we use the SB Galaxy model to constrain the age uncertainty. 
Although this model covers a wider range of ages, it covers a smaller range of metallicity, going no 
lower than around —2. At this metallicity, there is a negligible change in distance with a change in 
age. 
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Fig. 4. — Dwarf targets in various open and globular clusters compared to shifted YREC isochrones 
of comparable age (listed in Table[2]). For each cluster, we extract all targets with logg >4.1 and 
calculate their distance using a set of different methods. We have shifted the isochrones according 
to the mean distance for each cluster from the YREC isochrone technique (in black). The short 
dashed lines shift the isochrone by the la dispersion of the distance measurements for each cluster. 
The gray lines show the photometric parallax relationship using the mean SSPP metallicity of the 
clusters (see ^ IB.Sp . The vertical long-dashed lines show the SEGUE G- and K-dwarf color cut 
converted to {g — i). Each plot is labeled with the name of the cluster and the distance modulus 
used to shift the isochrones. The cluster distances determined using all of our methods are listed 
in Table[2l 
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Fig. 5. — The change in each [Fe/H] bin of the MDF due to different uncertainties. We run a 
Monte Carlo analysis on the TRILEGAL simulation convolving the expected uncertainty in each 
of these parameters with a Gaussian to examine how the structure of the MDF changes. At the 
metal-poor end, the number of stars in each bin is quite small. Thus, the uncertainty from different 
factors is very large. Typically, uncertainties from our bootstrap analysis and random distance 
errors dominate. The random distance uncertainties shown here neglect the correlated [Fe/H]- 
distance error. The MDF uncertainties presented in the following figures and tables reflect the 
total uncertainty per bin, combining all of the shown errors in quadrature. 
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Fig. 6. — A comparison of the surface gravity values for the grid of synthetic spectra of dwarfs 
and subgiants to those determined by the SSPP. These synthetic spectra were degraded to be at 
a S/N of 25, where the expected SSPP uncertainty in log (7 is ±0.40 dex. The top panel shows 
the difference between the model and SSPP-calculated surface gravities with respect to that of 
the model. The gray open circles show the grid points where modeled stars with dwarf surface 
gravities have SSPP log 5 of less than 4.1. The dashed lines show a change in log 5 of and ±0.4 
dex. The SSPP tends to underestimate the surface gravity of these synthetic spectra. There is some 
structure in A log g at the high log g end; this is likely due to the continuum removal processes in 
the SSPP and is being examined for the SEGUE DR9 public data release. The middle figure shows 
the change in surface gravity with respect to model effective temperature. Dwarfs misidentified as 
subgiants tend to be at hotter temperatures. Finally, the bottom panel shows the distribution of 
A log 5 with respect to model [Fe/H]. Misidentified dwarfs tend to be the most metal-rich stars, 
although there are a few occurences at the metal-poor end, where weak lines make it difficult to 
measure atmospheric parameters. The most metal-rich stars tend to have surface gravities close to 
the cutoff; even small underestimates can force them from the sample. 
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Fig. 7. — Comparison of the true cumulative metallicity distribution from the TRILEGAL model 
to a distribution which includes the correlated uncertainties in [Fe/H] and distance. We simulate 
the G (top, black) and K (bottom, black) dwarfs in SEGUE using a TRILEGAL model. Using a 
Monte-Carlo analysis, we factor in the uncertainties in [Fe/H] and the resulting change in distance. 
This broadens the metallicity distribution a significant amount, as indicated by the gray d([Fe/H]) 
lines. We estimate the slopes of these distributions between a cumulative fraction of 0.25 and 0.75; 
they are listed on the figure. This change estimates the amount of broadening in our MDFs caused 
by errors in metallicity. 
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Fig. 8. — The SEGUE G- and K-dwarf sample overlap with other SEGUE target categories, biasing 
them in metallicity space. The solid black line shows the ro distribution of all photometric stars 
along the SEGUE lines of sight that meet the color and magnitude criteria of a K dwarf. On the left 
hand side, we show the distribution of all photometric stars that meet only the criteria of a K dwarf 
as the dotted line. We also present the tq distribution of all SEGUE spectroscopic observations that 
meet only the criteria of a K dwarf. In contrast, on the right hand side we show the photometric 
and spectroscopic sample of stars that meet the criteria of a K dwarf and a low- metallicity target. 
It has a significantly different distribution than the K-dwarf targets alone. 
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Fig. 9. — A comparison of the mean target-type weight with respect to [Fe/H] for G (sohd Hne) and 
K dwarfs (dashed hne). There are fewer SEGUE fibers devoted to K-dwarf targets along each hne 
of sight, making the sample more contaminated by stars assigned fibers for other stellar categories. 
The target-type weights account for this by scaling down the proportion of metal-poor stars and 
scaling the metal-rich end of the distribution. The variation in target-type weights for G-dwarf 
stars is much smaller, as it suffers less from target-selection biases. 
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Fig. 10. — The raw and weighted MDF for the G (left) and K (right) dwarfs. The raw spectroscopic 
sample is the dashed line. The sample adjusted for target-type, r-magnitude, and mass-function 
weights is the solid line. The error bars on the solid line reflect the error in each bin based on our 
total MDF uncertainty, combining our bootstrap analysis with other errors (^ 15. 2|) . Both samples 
are limited to their spectral-type distance ranges, as described in § 14.51 Note that the corrections 
associated with the target selection scheme boost the metal-rich end of the distribution. 
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Fig. 11. — A test of the target-type and r-magnitude weighting system on the TRILEGAL Galaxy 
model. The solid black line represents the distribution of a randomly selected subsample of G and K 
dwarfs representing the SEGUE spectroscopic sample, which we refer to as the "raw spectroscopic." 
The gray line is the distribution of every modeled star that meets the G- and K-dwarf photometric 
criteria. Finally, the dashed line is the distribution of the spectroscopic sample once we have 
weighted it in r-magnitude and target type. The top figure shows the numerical distribution of the 
three samples; the good match between the weighted spectroscopic and photometric distributions 
indicates that our technique accurately corrects for SEGUE target selection biases. The bottom 
figure shows the fractional distribution of the three samples. The "raw" spectroscopic sample is 
biased towards metal-poor stars, but our weighting scheme corrects for this. 
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Fig. 12. — A test of the mass-function weights on the SEGUE hnes of sight modeled with TRILE- 
GAL. The top figure compares the effect of a color cut on the metallicity distribution to a mass 
cut. The color cut extracts all stars within the magnitude and {g — r) range of the SEGUE G and 
K dwarfs, shown in black. The mass cut selects all stars within the magnitude range with masses 
between 0.5 and 0.6 Mq, shown in gray. Each metallicity bin of the two samples probes the same 
volume of space. The two distributions are quite similar, although there are some discrepancies at 
the metal-poor end. The bottom row tests our weighting scheme. The gray distribution remains 
the same as in the top plot. The additional distributions represent a randomly chosen G/K dwarf 
subsample to represent SEGUE spectroscopy (black) and this distribution weighted in target type, 
r-magnitude, and mass function (dashed). The weighted distribution scales the raw spectroscopic 
down at the metal-poor end and boosts the metal-rich end, such that it simulates isolating a uniform 
portion of the mass function. 
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Fig. 13. — Comparison of the G- (black) and K-dwarf (gray) MDFs. Both samples are limited to 
distances between 1.59 and 1.84 kpc. The error bars are based upon our total MDF uncertainty, 
combining a bootstrap analysis over 500 iterations with other sample uncertainties (§ 15. 2p . The two 
samples cover a similar metallicity range and peak at approximately the same value. At low \Z\, 
K dwarfs have more metal-rich stars, due to evolutionary effect. At high \Z\, they also have more 
metal-poor stars. This results in a broader distribution than observed for the G dwarfs over all \Z\. 
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Fig. 14. — The spatial positions of our G (black) and K (red) dwarfs. 



Copenhagen Survey (GCS) sample of stars is shown in blue (iNordstrom et al. 



2008; 



Casagrande et al. 



2011 



The Geneva- 
Haywoodl 



2004 



The G dwarfs are limited to distances between 1.59 and 2.29 kpc, 
whereas K dwarfs fall between 1.18 and 1.84 kpc, based on the magnitude limits of SEGUE. This re- 
duces our stellar sample to 5,407 K- and 7,834 G-dwarfs. Our sample probes the disk both towards 
and away from the Galactic center and above and below the Galactic plane, with the line-of-sight 
structure of SEGUE clearly visible. The GCS sample is limited to the solar neighborhood, with 
distances typically less than 200 pc. 



PI 




Fig. 15. — The distribution of G (bottom) and K dwarfs (top) in [Fe/H] with respect to \Z\. Stars 
are sorted into bins of [Fe/H] and \Z\. The number of stars in each bin is then adjusted using the 
target-type, r-magnitude, and mass-function weights. The logarithm of the weighted number of 
targets is indicated by color, as labeled on the right side. 
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Fig. 16. — The weighted MDF over the overlap distance range for the G (black) and K (gray) 
dwarfs at a range of \Z\, where |^| is the current distance from the plane in kpc. The uncertainty 
in each bin was determined from our bootstrap analysis over 500 iterations combined with additional 
measured uncertainties (§ 15. 2p . These distributions are listed in Table[6] and [71 
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Fig. 17. — The cumulative weighted MDF over the overlap distance range for the G (black) and 
K (gray) dwarfs at a range of |Z|, where \Z\ is the current distance from the plane in kpc. We 
combine the error bars for each individual metallicity bin (from our bootstrap and Monte Carlo 
analysis) together in quadrature to determine the uncertainty for each cumulative bin. Above \Z\ 
of 0.5 kpc, the K dwarfs show a more prominent metal-poor tail than the G dwarfs. Note that the 
fractions of the G- and K-dwarf sample below [Fe/H] of —1.0 are consistent until \Z\ >1.5 kpc. 
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Fig. 18. — Comparison of the G- (left) and K-dwarf (right) sa mple in bins of \Z\ with a simple closed 
box model (dotted line) assuming instantaneous recycling (jSchmidtlll963l ). The total metallicity 
distribution for each spectral type over its associated distance range is shown in the top row. The 
size of the discrepancy between the observed distribution and the model at the metal-poor end is 
consistent for G and K dwarfs. The G- and K-dwarf problem is evident throughout the Milky Way 
disk. 
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Fig. 19. — Compa rison of the SEGUE G- (left) and K-dwarf (right) MDFs with the distribution of 



Katz et al.l (|201ll ). The S EGUE samples ar e represented by the solid line, they cover the spectral- 



type d istance ranges. Th e lKatz et alJ (|201ll ) distributions are plotted as the dashed lines. The MDF 
of the iKatz et al.l (120111 ) sample is typically more metal-rich than that of the G and K dwarfs. 
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Fig. 2 0. — Comparison o f the SEGUE G-dwarf MDF over the spectral-type distance range with 
that of iLee et al.l (j2011bl ). Their work is based on the SEGUE G-dwarf sample and di vides their 



stars by [a/Fe]. Our SEGUE G dwarfs are plotted in black. The blue dashed line is the lLee et al, 



(|2011bl ) a-poor sample, as a fraction of their entire sample of G dwar fs. The re d dashed line is their 
a-enhanced distribution and normalized in the same manner. The iLee et al.l (|2011bl ) sample has 
not been corrected for observational biases; thus, the two distributions are skewed to be more metal- 
poor. At small \Z\^ our corrected distributions appears to sample both the a components. The 
MDFs of our corrected distributions shifts more metal-poor with increasing |Z|, better matching 
the a-enhanced sample. 
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Fig. 21. — Comparison of the distribution of G-dwarf stars in the SEGUE sample (sohd black) 
to that in the SB Galaxy models (dashed black) in the left column. Both sample and model are 
limited to distances between 1.59 and 2.29 kpc. The SB Galaxy model is significantly more metal- 
rich than the SEGUE G-dwarf sample. This discrepancy is also seen for the K-dwarf sample. On 
the right side, we com pare the G - dwarf distribution in the model (dashed black) to the a-separated 
G-dwarf sample from lLee et al.l (|2011bl ). The sample with enhanced [a/Fe] is shown in red; that 
with low [a/Fe] is plotted in blue. For \Z\ less than 1 kpc, the SB model is consistent with the 
Qf-poor distribution. Above 1 kpc, the model is more metal-rich than this component. The general 
agreement with solely the a-poor portion suggests that the SB models currently only recreate the 
chemistry of the this portion of the Galaxy. 
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Fig. 22. — The sam e as Fig urel2 H except now for the TRILEGAL models, rather than those of 
Schonrich BinnevI (|2009al lb[) . Note that TRILEGAL simulates total metallicity, [M/H], rather 



than [Fe/H]. We assume a simple linear relationship between [Fe/H] and [a/Fe] to adjust the 
model from [M/H] to [Fe/H]. This model predicts more metal -rich MDF s at all heights than our 
observed sample. Comparison with the a-separated samples of lLee et al.l (j2011bl ) suggests that at 
high |Z|, TRILEGAL does not adequately model the a-rich component. However, we can improve 
the agreement by adjusting the assumed [a/Fe]-[Fe/H] relationship. The comparison of SEGUE K 
dwarfs with the TRILEGAL model shows similar behavior. 
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Fig. 23. — The surface gravities calculat ed using different methods in the SSPP for our sample 
of G and K dwarfs (see Lee et al. 2008J jbl for a description of methods). The original adopted 
log 5 distribution is in the top left. All of the distributions have a tail extending to logg' values 
indicative of giants. However, this tail is significantly more pronounced for the CaI2 and MgH 
methods. Additionally, the width of the peaks for these two methods are larger than that of the 
others, making their estimate of log 5 generally more uncertain. The WBG method is, in general, 
unable to calculate a surface gravity for our targets. When the CaI2, MgH, and WBG techniques 
are removed, the low log 5 tail is diminished, and the peak of the distribution is narrower. This 
corrected log (7 is shown at the top right. 



- 69 - 




Fig. 24. — A comparison of the atmospheric parameters determined by the optimized version of 
the SSPP with those from DR8. The left column directly compares the values, with the optimized 
version on the x axis and the DR8 on the y. The red line shows a one-to-one correlation. The right 
column shows the distribution of the difference between the two values for each parameter, with 
the mean and standard deviation of the difference noted in the top left corner. The top row is for 
[Fe/H] and bottom for logg. 
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Fig. 25. — The metallicities calculated using different methods in the SSPP for our G- and K-dwarf 
sample. The original adopted [Fe/H] value is on the top left, with the adjusted distribution on the 
top right. For the different methods used in our optimized version of the SSPP, we plot the original 
distribution for the sample in black and the modified version in red. Note that for many methods, 
such as NGSl and CallKl, the modified distribution is much smoother. We have also eliminated a 
number of methods as they are either not applicable for our sample or do not cover the appropriate 
metallicity range. 
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Table 1. Cluster Metallicities determined with the Optimized SSPP 







[Fe/H] 




Cluster 


Literature 


All 


(.9 - r) 


M92 


-2.25 


-2.23±0.22 


-2.32±0.09 


M3 


-1.55 


-1.55±0.18 


-1.58±0.16 


M71 


-0.79 


-0.78±0.12 


-0.77±0.13 


NGC 2420 


-0.20^^ 


-0.28±0.13 


-0.26±0.13 


NGC 2158 


-0.26 


-0.26±0.10 


-0.30±0.09 


M67 


O.OS'' 


0.06±0.07 


0.03±0.10 


NGC 6791 


0.31 


0.40±0.12 


0.36±0.09 



Jacobson et al.l (2011) 



iRandich at al.l (|2006l ) 



Note. — The mean and a [Fe/H] values of clus- 
ter members based upon the optimiz ed version of the 
SSPP . All literature values are from ISmolinski et al.l 
POIII ) unless otherwise noted. The calculated metal- 
licities are similar to the literature values demonstrat- 
ing that this SSPP version accurately determines stel- 
lar parameters. The "((7 — r)" column isolates mem- 
bers with colors in the appropriate range for G and 
K dwarfs. These values are not significantly different 
than those for the whole cluster sample. 



Table 2. Cluster Distances 



Cluster 


[Fc/H] 


E(B-V) 


Assumed Age 


Literature 


Ivezic et al. (2008) Photometric tt 


YREC 








Gyr 


kpc 


kpc 


kpc 


M13 


-1.54 


0.017 


14 


7.7 


7.71+0.78 


8.01+1.01 


M67 


+0.02 


0.032 


4 


0.91 


0.93+0.10 


0.87+0.12 


NGC 2420 


-0.44 


0.041 


2 


3.09 


2.73+0.26 


2.39+0.19 


NGC 6791 


+0.30 


0.117 


4 


4.10 


3.83+0.69 


3.47+0.52 



Not e. — The par ameters a n d dista nces determined for the test clusters. The E(B-V) is based on Schlegel et al 



(jl998l ) , as listed in iLee et al 



(|2008bl) . The listed metallicity and literature distances arc from iHarrisI (119961 ) for 
the three globular clusters and WEBDA for the two open clusters. Note that these metallicity values are slightly 
different than those found using the SSPP listed in Tabled the [Fe/H] hsted in this table is estimated using the 
optimized version of the SSPP (§ lA| . 



Table 3. Total Metallicity Distribution Functions of G and K dwarfs over Different Distance Ranges 



G dwarfs K dwarfs 

[Fc/H] 1.59<d<2.29 kpc 1.59<d<1.84 kpc 1.18<d<1.84 kpc 1.59<d<1.84 kpc 





Raw 


Weighted 


Raw 


Weij 


^hted 


Raw 


Weij 


;hted 


Raw 


Weighted 


-3.3 


-0 





4-0 
-0 


-0 





+0 

-0 


1+1 

-1 


7 


+6 
-7 


-0 


» ^i; 

-0 


-3.1 




1 


+1 

-1 


-0 





+0 

-0 


« +2 
6 o 

-3 


25 


+49 
-13 


-2 


18 

-11 


-2.9 


-2 


15 


+15 
-10 


-0 





+0 
-0 


-2 


22 


+17 
-10 






-2.7 


-3 


26 


+10 
-13 


-2 


10 


+6 
-7 


16 A 

-4 


58 


+83 
-21 




1« ""In 

-10 


-2.5 


-3 


9 


+7 

-3 


-2 


4 


+3 

-3 


-5 


57 


+22 
-17 


-2 


18 


-2.3 


-5 


38 


+24 
-11 


-4 


15 


+7 
-6 


-6 


507 


+206 
-149 


20+^ 

-4 


-136 


-2.1 


+14 
-6 


190 


+80 
-53 


^°^!: 

-4 


33 


+33 
-13 


-9 


673 


+224 
-150 


-7 


-167 


-1.9 


-10 


110 


+28 

1 8 


25+' 
-4 


24 


+17 

-6 


-12 


1151 


+780 
-238 


fir; +11 

65 n 

-9 


298 «f 

-52 


-1.7 


+21 


388 


+134 
-126 


+ 12 
^0-8 


95 


+47 
-26 


+27 
'"-13 


1492 


+670 

-440 


+ 14 
62 g 


+413 
-455 


-1.5 




425 


+151 

-68 


fir, +15 
65 _8 


165 


+70 

-47 




1880 


+346 
-306 






-1.3 




611 


+169 
-113 




238 


+120 

-94 




1637 


+472 
-248 






-1.1 




1324 


+326 
-241 




390 


+135 
-90 




3074 


+884 
-420 






-0.9 




2083 


+404 
-163 


-22 


622 


+218 
-82 


-!2^ 


4253 


+545 
-407 


-19 


20.«- 


-0.7 


1700+"'^ 


1.072x10 


4+1027 


645+^^ 


3618 


+540 


908+^1 


1.309x10^+1^^1 


333+2^ 


5941+™'' 



Table 3 — Continued 



G dwarfs K dwarfs 



[Fe/H] 1.59<d<2.29 kpc 1.59<d<1.84 kpc 1.18<d<1.84 kpc 1.59<d<1.84 kpc 





R,cLW 


Weighted 




V V ClgllUOJ. 




crVifpr] 
V V CitlllJV^U. 


R,cLW 


V V d^llLCU. 




-56 




-707 


-30 


-411 


-32 


-848 


-20 


-888 


-0.5 


+77 
1550_^2 


1.62x10^ 


+1195 
-691 












+1228 


-0.3 


+24 
1045_^3 


1.146x10^ 


+655 
-776 














-0.1 


+27 
^20-28 


6007 


+447 
-378 














0.1 


+ 19 


1258 


+115 
-186 








-Z 






0.3 




168 


+30 

-28 


45 _^ 






1» 






0.5 







+0 

-0 















Note. — The number of G and K dwarfs in each metallicity bin using different distance hmits. The Raw numbers 
are the original spectroscopic sample; the Weighted are once we have applied our corrections that accounting for 
SEGUE target selection effects. The listed errors reflect the total uncertainties in each bin determined by combining 
a bootstrap analysis over 500 iterations with Monte Carlos over various properties (§[5?2|). 
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Table 4. The Metallicity Distribution Functions of G dwarfs over Different Ranges of |Z| for 

1.59<d<2.29 kpc. 



[Fe/H] 


0.0< 


Z| <0.5 kpc 


0.5< 


Z\ <1.0 kpc 


1.0< |Z| <1.5 kpc 


1.5< \Z\ <2.5 kpc 




Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


-3.3 


0+° 
-0 


+° 
-0 





-1-0 
-0 


+° 
-0 








° -S 


-3.1 


0+" 
-0 


+° 
-0 


1 


+1 

-1 


1 +1 

-1 








° -2 


-2.9 


0+° 
-0 


+° 
-0 


1 


-1-1 
-1 


12 +13 
-13 








3 ^2 


-2.7 


0+° 
-0 


+° 
-0 


2 


+1 
-1 


13 +12 
-9 


4-2 


9 -J 


2^? 


4 -3 


-2.5 


1+1 

-1 


+° 
-0 


2 


+2 
-2 


2 +2 
-2 


3^2 


3 ^2 


3 3 




-2.3 


1+1 

-1 


3 +^ 
-2 


3 


+3 
-2 


10 +24 
-6 




- 


12!: 


« 


-2.1 


1+1 

-1 


+° 
-0 


7 


+4 
-3 


103+'^ 
-52 


20 


- -22 


31^9 


29 -10 


-1.9 


0+° 
-0 


+° 
-0 


8 


+4 
-3 


10 +"5 

-4 


3^-9 


- -?9 




3^ -n 


-1.7 


3+2 
-2 


146+1^° 
-135 


31 


+8 
-8 


+39 

88 

-25 


--lo 




52+- 




-1.5 


1+3 
-1 


58 +1=^^ 
-50 


33 


+12 
-5 


201+139 
-48 


— O 


- ^8 


— 1 1 


- "16 


-1.3 


2+2 
-2 


121+1^^ 
-121 


32 


+1 1 

-6 


170+130 

-50 










-1.1 


— o 




87 


+28 
-13 




-'IS 








-0.9 






121 


+74 
-26 












-0.7 






286 


+30 
-25 












-0.5 






358 


+32 
-16 




+44 
^6°-42 




497+'= 

-24 




-0.3 




'-IS? 


350 


+15 
-24 




421+15 
-32 








-0.1 






318 


+17 
-21 












0.1 


< 




169 


+12 
-18 












0.3 






52 


+7 
-10 












0.5 






1 


+1 
-1 













Note. — The G-dwarf raw and weighted metallicity distribution function for different ranges of \Z\. This 
is for the G dwarf sample limited to distances between 1.59 and 2.29 kpc. 
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Table 5. The Metallicity Distribution Functions of K dwarfs over Different ranges of |Z| for 

1.18<d<1.84 kpc. 



[Fe/H] 


0.0< 


Z| <0.5 kpc 


0.5< |Z| <1.0 kpc 


1.0< \Z\ <1.5 kpc 


1.5< \Z\ <2.5 kpc 




Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


-3.3 


0+° 
-0 


+° 
-0 




7 


-1-6 
-7 




« -S 


°-0 




-3.1 


0+° 

-0 


+° 
-0 


it; 


6 


-1-46 
-6 








--8 


-2.9 


1 +1 

-1 


10 

-10 







-1-0 
-0 




. !: 




^-5 


-2.7 


1 +1 

-1 


21 

-15 


3-2 


21 


-1-79 
-11 










-2.5 


1 +1 

-1 


4 +^ 

-3 




27 


+18 
-14 










-2.3 


-1 


21 +6^ 

-14 


--5 


406 


-1-186 
-168 
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51 


-1.7 
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-3 


171+^^ 
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-1-643 
-448 








-^26 


-1.5 


12+^ 
-4 


337+^^^ 
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1019 


-1-292 
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-^"13 




22+5 




-1.3 
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-4 
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-^^^ 


1026 
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— ooo 




1078 


-1-641 
-193 






33+;° 




-0.9 
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Note. — The K-dwarf raw and weighted metallicity distribution function for different ranges of \Z\. This is 
for the distances between 1.18 and 1.84 kpc. 
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Table 6. The Metallicity Distribution Functions of G dwarfs over Different Ranges of \Z\ for 

1.59<d<1.84 kpc. 



[Fe/H] 


0.0< 


Z| <0.5 kpc 


0.5< \Z\ <1.0 kpc 


1.0< \Z\ <1.5 kpc 


1.5< \Z\ <2.5 kpc 




Raw 


Weighted 


Raw 


Weighted 


Raw 


Wei, 


ghted 


Raw 


Weighted 


-3.3 


0+° 
-0 


+° 
-0 











+0 
—0 






-3.1 


0+° 
-0 


+° 
-0 











+0 






-2.9 


0+° 
-0 


+° 

-0 




° -0 







+0 
—0 






-2.7 


0+° 

-0 


+° 

-0 








6 


+6 
_5 






-2.5 


1+1 

-1 


+° 
-0 




1 -? 




3 


+3 
_3 


1-? 


°-0 


-2.3 


1+1 

-1 


3 +^ 

-3 




3 -2 


— o 


9 


+6 

— o 




°"o 


-2.1 


1+1 

-1 


+° 
-0 


— u 


° 

— u 


1°1 


27 


+n 


'-r 


^1 

— o 


-1.9 


0+° 

-0 


+° 

-0 


«-3 


« -3 




14 


-hll 

—6 




^-3 


-1.7 


2+2 

-2 


16 +3^ 
-15 




- -2? 




18 


4-18 
_7 


11 -I 


^^3 


-1.5 


0+° 
-0 


+° 
-0 


— u 


i-'ei 

— uu 




24 


+17 
— u 


-^10 


111 


-1.3 


2+1 
-2 


121+132 
-98 




" -18 


+10 
^^-12 


57 


-1-23 
-20 






-1.1 




79 +11^ 
inn 






97+^° 
-32 


86 


+57 
-30 






-0.9 












242 


-1-431 
-248 






-0.7 


< 




-1-24 






814 


+112 
-172 






-0.5 


< 








-IS 


133S 


,-H155 
■-113 






-0.3 




-IS 








652 


+61 
-65 






-0.1 












297 


+28 
-40 






0.1 




-65 








52 


+11 
-12 






0.3 












12 


+5 

-4 






0.5 




«!2 











+0 
-0 







Note. — The G-dwarf raw and weighted motaUicity distribution function for different ranges of \Z\. This 
is for the distances between 1.59 and 1.84 kpc. 
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Table 7. The Metallicity Distribution Functions of K dwarfs over Different ranges of |Z| for 

1.59<d<1.84 kpc. 



[Fe/H] 


0.0< 


Z\ <0.5 kpc 


0.5< |Z| <1.0 kpc 


1.0< |Z| <1.5 kpc 


1.5< |Z| <2.5 kpc 




Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


Raw 


Weighted 


-3.3 


0+° 
-0 


+° 
-0 














-3.1 


-0 


+° 
-0 




, .49 


1-? 


1 -I 




--8 


-2.9 


-0 


+° 
-0 














-2.7 


-0 


+° 
-0 


2-2 


- -13 


1-1 


^ -f-19 






-2.5 


-0 


+° 
-0 




8+- 


^-3 


11 i 






-2.3 


1+1 

-1 


3 +=^^ 
-3 




267+1^1 
—148 




— yj 


*-2 




-2.1 


-2 


104+^3 
-82 


11-4 


279 +1^° 
-134 




- -26 






-1.9 


0+° 

-0 


+° 

-0 


18+^ 


142+398 
—49 




l°^-25 


--5 




-1.7 


1+1 

-1 


23 +=^3 
-19 


-2 


--:?3 




--25 






-1.5 


2+2 
-2 


65 +11^ 
-44 


— y 


-«"f76 

— 1 / u 




i-":8 




74+23 

9Q 


-1.3 


3+2 
-2 


88 +1*^^ 
-65 










22+8 
-5 




-1.1 




^1*^ 327 




-IS 










-0.9 








™^?^^ 










-0.7 


















-0.5 


















-0.3 














< 


-1-47 
'"■-42 


-0.1 


















0.1 










< 


314+93 
31* -108 






0.3 














< 




0.5 












«!2 







Note. — The K-dwarf raw and weighted metallicity distribution function for different ranges of \Z\. 
This is for the distances between 1.59 and 1.84 kpc. 



